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Beller,  Denis  Eugene,  Capt ,  USAF. 


Conceptual  Design  and  Neutronics  Analyses  of  a  Fusion 
Reactor  Blanket  Simulation  Facility. 


183  pages,  Ph.D.,  Purdue  University,  1986.  Major 
Professor:  K.O.  Ott. 

A  new  conceptual  design  of  a  fusion  reactor  blanket 
simulation  facility  has  been  developed.  This  design 
follows  the  principles  that  have  been  successfully 
employed  in  the  Purdue  Fast  Breeder  Blanket  Facility 
(FBBF),  because  experiments  conducted  in  it  have 
resulted  in  the  discovery  of  deficiencies  in 
neutronics  prediction  methods.  With  this  design, 
discrepancies  between  calculation  and  experimental 
data  can  be  fully  attributed  to  calculation  methods 
because  design  deficiencies  which  could  affect  results 
are  insignificant.  Inelastic  scattering  cross 
sections  are  identified  as  a  major  source  of  these 
discrepancies.  The  conceptual  design  of  this  FBBF- 
analog,  the  fusion  reactor  blanket  facility  (FRBF),  is 
presented.  Essential  features  are  a  cylindrical 
geometry  and  a  distributed,  cosine-shaped  line  source 


of  14  MeV  neutrons.  This  source  can  be  created  by 
sweeping  a  deuteron  beam  over  an  elongated  titanium- 
tritide  target.  To  demonstrate  that  the  design  of  the 
FRBF  will  not  contribute  significant  deviations  in 
experimental  results,  neutronics  analyses  were 
performed:  results  of  comparisons  of  2— dimensional  to 
1-dimensional  predictions  are  reported  for  two  blanket 
compositions.  Expected  deviations  from  1-D 
predictions  which  are  due  to  source  anisotropy  and 
blanket  asymmetry  are  minimal.  Then  1-D  calculations 
can  be  performed  in  fine  detail  to  produce  course- 
group  constants  for  2-D  predictions.  When  these  2-D 
neutronics  studies  include  the  asymmetry  and 
anisotropy,  differences  between  prediction  and 
experimental  results  caused  by  the  design  will  be  less 
than  1SS.  Thus,  the  design  of  the  FRBF  allows  simple 
and  straightforward  interpretation  of  experimental 
results,  without  a  need  for  course  3-D  calculations. 
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Belter,  Denis  Eugene.  Ph.D.,  Purdue  University,  May  1986. 
Conceptual  Design  and  Neutronics  Analyses  of  a  Fusion 
Reactor  Blanket  Simulation  Facility.  Major  Professor:  K.O. 
Ott. 

A  new  conceptual  design  of  a  fusion  reactor  blanket 
simulation  facility  has  been  developed.  This  design  follows 
the  principles  that  have  been  successfully  employed  in  the 
Purdue  Fast  Breeder  Blanket  Facility  (FBBF) ,  because 
experiments  conducted  in  it  have  resulted  in  the  discovery 
of  deficiencies  in  neutronics  prediction  methods.  With  this 
design,  discrepancies  between  calculation  and  experimental 
data  can  be  fully  attributed  to  calculation  methods  because 
design  deficiencies  which  could  affect  results  are 
insignificant.  Inelastic  scattering  cross  sections  are 
identified  as  a  major  source  of  these  discrepancies.  The 
conceptual  design  of  this  FBBF— analog,  the  fusion  reactor 
blanket  facility  (FRBF),  is  presented.  Essential  features 
are  a  cylindrical  geometry  and  a  distributed,  cosine-shaped 
line  source  of  14  MeV  neutrons.  This  source  can  be  created 
by  sweeping  a  deuteron  beam  over  an  elongated  titanium- 
tritide  target.  To  demonstrate  that  the  design  of  the  FRBF 
will  not  contribute  significant  deviations  in  experimental 
results,  neutronics  ah^lyses  were  performed:  results  of 
comparisons  of  2— dimens i/onal  to  1— dimensional  predictions 
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are  reported  for  two  blanket  compositions.  Expected 
deviations  from  1-D  predictions  which  are  due  to  source 
anisotropy  and  blanket  asymmetry  are  minimal.  Then  1-D 
calculations  can  be  performed  in  fine  detail  to  produce 
course-group  constants  for  2-D  predictions.  When  these  2-D 
neutronics  studies  include  the  asymmetry  and  anisotropy, 
differences  between  prediction  and  experimental  results 
caused  by  the  design  will  be  less  than  IX.  Thus,  the  design 
of  the  FRBF  allows  simple  and  straightforward  interpretation 
of  experimental  results,  without  a  need  for  course  3-D 


calculations. 


1 


I.  INTRODUCTION 

Past  work  in  fusion  reactor  blanket  neutronics  indicates 
a  stronf  need  for  an  inexpensive,  flexible,  properly 
designed  experimental  facility.  This  is  pointed  out  by  the 
FINESSE  study1 ,  a  recent  report  on  future  requirements  and 
plans  for  fusion  reactor  blanket  and  materials  testing.  It 
expresses  the  importance  of  accurately  predicting,  within  at 
least  5%,  neutron  transport  and  tritium  breeding  rates.  The 
reasons  why  accurate  prediction  of  the  transport  of  neutrons 
through  the  first  wall,  blanket,  and  shielding  is  of  such 
great  importance  in  future  designs  of  fusion  reactors  is 
because  in  DT-fueled  reactors  1)  the  neutrons  carry  80 
percent  of  the  reaction  energy  away  from  the  plasma,  2)  they 
cause  damage  to  and  activition  of  reactor  materials,  and  3) 
they  must  provide  new  fuel  through  breeding  reactions.  A 
major  source  of  error  in  these  predictions  will  be  from 
current  uncertainties  in  inelastic  scattering  matrices  in 
cross  section  libraries.  Approximately  75-90X  of  the  14  MeV 
of  DT-neutron  energy  is  deposited  as  a  result  of  inelastic 
scattering  reactions.  (see  section  V.3,  Cross  Sections) 
Thus,  errors  in  inelastic  data  can  cause  significant  errors 
in  predicting  both  the  location  and  amount  of  neutron- 


induced  energy  (gamma  and  other  products)  deposition  in 
proposed  fusion  reactor  designs. 

However,  the  few  experiments  which  were  done  to  compare 
to  transport  predictions  in  simulated  fusion  reactor 
blankets  have  not  produced  data  which  validates  the 
neutronics  methods  or  the  data  libraries.  The  facilities 
which  had  a  "clean"  design  were  inflexible;  thus  they  could 
not  produce  multiple  experiments  to  compare  to  a  given 
neutronics  method.  More  experiments  would  require 
construction  of  a  new  system,  at  great  expense.  Generally, 
the  results  of  these  individual  experiments  deviated  from 
predictions  by  20  to  40*. 2-8  The  results  of  experiments 
done  in  flexible  facilities  deviated  up  to  50*  from 
calculations.  In  only  a  few,  one-of-a-kind  experiments  were 
results  produced  which  compared  favorably  (about  10* 
deviations)  to  one  of  a  number  of  prediction  methods. 

Transport  codes  are  generally  used  in  fusion  reactor 
design  studies  to  evaluate  tritium  breeding,  heat 
deposition,  and  neutron  leakage  to  field  magnets;  and  design 
decisions  are  based  upon  these  evaluations.  Confidence  in 
the  results  of  these  studies  can  only  be  achieved  by 
improving  the  transport  codes  and  group  constants  used  to 
predict  neutron  behavior  and  by  experimental  validation  of 
these  codes.  Generally,  determination  of  the  causes  of 
error*  in  past  transport  predictions  has  been  very  difficult 
because  of  the  system  geometries  and  the  neutron  sources 
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used  in  the  experiments.  These  "deficiencies”  of  the  set¬ 
ups  clearly  affect  the  validation  of  codes  and  high-energy 
cross  sections.  Eliminating  the  interpretation  problems  due 
to  deficiencies  in  the  design  of  the  fusion  reactor  blanket 
simulation  facilities  is  the  subject  of  this  research. 

Qualitatively,  the  effects  of  the  deficiencies  in  these 
experimental  set-ups  should  be  considerably  less  than  the 
acceptable  deviations  between  eventual  experimental  results 
and  theoretical  predictions.  This  is  a  guideline  which  can 
be  used  during  evaluation  of  conceptual  designs  of  these 
facilities;  however,  this  guideline  has  not  been  met  for 
past  fusion  reactor  blanket  neutronics  evaluations. 

Fusion  reactor  blanket  experiments  are  conducted  to 
determine  the  space  and  energy  dependence  of  neutron  flux 
and  reaction  rates.  Because  of  the  computing  time  and 
expense  required  for  a  theoretical  three-dimensional  (3— D) 
description  of  flux  spectra,  attempts  are  made  to  reduce  the 
space  dependencies  (or  to  simplify  their  description  in  one 
or  two  dimensions).  If  an  appropriate  geometry  and  incident 
neutron  flux  is  used,  the  space  dependency  can  essentially 
be  reduced  to  one  dimension.  A  slightly  less  advantageous* 
but  still  very  useful  technique,  is  to  use  a  geometry  and  a 
source  that  allow  separation  of  the  flux  function  into  a 
spatial  dependency  in  one  dimension  and  a  space/energy 
dependency  in  the  other  dimension.  Then  the  lateral 
transport  can  be  approximated  analytically  while  the  space 
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and  enargy  variation  can  be  treated  in  one  dimension  using 
sophisticated  codes.  While  not  exact,  this  separability  is 
a  very  good  approximation,  and  the  resulting  differences  are 
less  than  the  deviations  being  investigated.  However, 
previous  fusion  reactor  blanket  experiments  have  not 
provided  the  geometry  and  source  necessary  to  permit  these 
simplified  computations  without  introducing  corrections 
which  were  greater  than  the  effects  being  investigated. 

Past  simulated  fusion  reactor  blanket  experiments  have 
been  designed  with  parallelepiped  (slab),  spherical,  or 
cylindrical  geometries. These  set-ups  have  been 
deficient  in  various  ways: 

1)  An  infinite  slab  with  a  plane-wave  neutron  source  is 
the  ideal  geometry  for  1-D  neutronics  analysis.  However, 
the  thick  parallelepiped  slabs  used  in  experiments  do  not 
approximate  infinite  slabs,  and  the  neutron  flux  from  a 
smal 1  target  is  not  a  plane— wave  source  (nor  does  it  impose 
the  proper  boundary  condition  for  a  non-infinite  slab). 
Lateral  leakage  dominates  the  transport  in  these  slabs;  it 
is  accounted  for  by  including  a  buckling  term  in  the 
neutronics  codes  or  by  approximating  the  lateral  transport 
in  radial  geometry.  Errors  in  buckling  corrections  can 
overshadow  other  errors  in  the  neutronics  calculations,  and 
a  square  slab  cannot  be  modelled  in  radial  geometry. 

2)  A  system  with  spherical  geometry  and  a  centered 
isotropic  point  source  can  also  be  modelled  in  1—0  geometry. 
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However,  a  DT— neutron  generator  which  produces  ~14  MeV 
neutrons  was  used  for  these  experiments.  These  neutron 
sources  are  anisotropic  in  energy  and  yield,  which 
eliminates  the  advantage  of  the  spherical  geometry  (although 
the  system  does  retain  azimuthal  symmetry).  Two  other  major 
disadvantages  of  this  geometry  are  its  total  lack  of 
flexibility  to  modify  the  blanket  composition  and  the 
requirement  to  bore  asymmetric  holes  for  experiments  and 
measurements  inside  the  blanket. 

3)  Another  geometry  which  allows  1—0  analysis  is  an 
infinite  cylinder  with  an  axial  line  source.  Alternately,  a 
finite-length  cylinder  with  a  cosine-shaped  source 
distribution  as  a  boundary  condition  in  the  center  axis  can 
be  analyzed  with  1-D  calculations  in  good  approximation. 
However,  the  past  experiments  which  used  cylindrical 
geometry  used  anisotropic  point  sources  at  the  center  of  the 
■rid  plawe.  Again  the  lateral  flux  adjustment  and  leakage- 
differences  in  this  geometry  strongly  affect  the  transport, 
which  makes  this  set-up  inappropriate  for  validation  of 
codes.  These  deficiencies  in  the  design  of  the  experiments 
and  other  problems  resulted,  almost  without  exception,  in 
overestimation  of  tritium  breeding  and  incorrect  spectral 
predictions,  with  discrepancies  from  5  percent  to  50 
percent . 

These  fusion  reactor  blanket  simulation  facilities  can  be 
contrasted  with  the  Purdue  Fast  Breeder  Blank  t  Facility 


Many 


(FBBF) ,  which  allows  very  clean  experiments.*1 
different  kinds  of  measurements  have  been  made  in  the  FBBF, 
and  all  consistently  show  the  same  trends.  Although  the 
neutronics  problems  are  milder,  due  to  lower-energy  fission 
neutrons,  discrepancies  have  also  been  found  with  FBBF 
experiments.  Unlike  the  fusion  reactor  blanket  experiments, 
bhese  differences  can  be  fully  attributed  to  deficiencies  in 
the  computational  methods  and  cross  sections,  and  not  to 
deficiencies  in  the  configuration  of  the  experimental 
facility.  This  is  because  the  cylindrical  facility  has  a 
distributed  neutron  source  in  the  axis  with  a  cosine— shaped 
flux  distribution.  Axial  flux  curvature  which  develops  in 
the  cylindrical  blanket  is  already  imposed  by  the  source; 
thus  little  adjustment  is  required.  Therefore,  radial 
transport  can  be  studied  while  errors  due  to  axial  leakage 
are  minimized.  Two  other  features  contribute  to  the 
improved  treatment  of  the  neutron  flux:  the  FBBF  has 
complete  azimuthal  symmetry,  and  the  ratio  of  height  to 
radius  is  great  enough  that  the  flux  curvature  is  minimal. 
Thus,  there  is  no  azimuthal  dependence  and  almost  no  axial 
dependence  of  the  flux  at  the  mid-plane  of  the  blanket. 

These  features  have  allowed  researchers  working  with  the 
FBBF  to  determine  the  magnitudes  of  errors  in  neutron 
transport  predictions  with  unprecedented  accuracy,  and  to 
determine  indications  of  their  causes.  It  has  also 
permitted  them  to  discover  differences  between  measurement 
techniques  and  to  determine  the  sources  of  errors  in  some  of 


them.  The  sam*  should  be  expected  of  a  fusion  reactor 
blanket  simulation  facility  based  on  the  same  principles. 

The  research  reported  in  this  paper  examines  the 
conceptual  design  of  a  cylindrical  facility  similar  to  the 
FBBF:  the  Fusion  Reactor  Blanket  Facility.  This  special 
facility  will  be  referred  to  henceforth  in  this  report  as 
the  FRBF.  The  FRBF  will  have  the  same  general  purpose  as 

the  FBBF,  ” . to  allow  a  direct  comparison  of  measurements 

with  design  code  calculations  in  a  realistic  geometrical 
structure.’81  The  FRBF  will  be  cylindrical,  to  simulate  a 
section  of  the  blanket  of  a  toroidal  or  mirror  reactor  and 
to  take  advantage  of  simplified  geometry;  it  will  be 
flexible,  to  allow  simulation  of  various  blanket 
configurations  and  materials;  it  will  have  a  cosine-shaped 
axially  distributed  neutron  source,  to  greatly  increase  the 
accuracy  of  a  buckling  correction;  it  will  take  advantage  of 
previously  developed  and  thoroughly  tested  neutronics  and 
reaction  rate  measurement  techniques  used  in  the  FBBF;  and 
it  will  take  advantage  of  a  faculty  and  staff  experienced  in 
the  measurement  of  high-energy  neutrons,  theoretical 
transport  predictions,  and  experiments  in  breeder  blankets. 
Studies  conducted  in  the  FRBF  will  greatly  enhance  the 
expanding  body  of  methodologies  for  fusion  reactor  blanket 
design,  and  will  add  confidence  to  studies  leading  to  future 
fusion  reactor  design  decisions. 
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The  remainder  of  this  paper  befins  with  a  review  of 
fusion  reactor  blanket  research  in  Chapter  II,  specifically 
a  review  of  past  blanket  experiments  and  comparisons  of 
experimental  results  to  calculations.  An  overview  of  the 
FBBF  is  presented  in  Chapter  III,  including  a  brief 
description  of  the  facility  and  results  of  FBBF  experiments. 
Then,  the  design  of  the  FRBF  itself  is  described  in  Chapter 
IV,  including  discussions  of  the  basis  for  selecting  the 
blankets  to  be  studied,  of  the  anisotropy  of  the  neutron 
source,  and  of  the  requirements  for  analysis  of  the  design. 
Chapter  V  includes  a  discussion  of  the  transport  code  (DORT) 
used  to  examine  neutron  transport,  tritium  breeding,  and 
leakage;  it  includes  a  discussion  of  the  modelling  of  the 
blankets  (1-D,  symmetric  2-D,  and  asymmetric  2-D) ;  and  it 
includes  a  discussion  of  the  treatment  of  neutron  sources 
(isotropic  point,  anisotropic  boundary,  and  distributed 
plasma  sources).  The  results  of  the  transport  computations 
are  presented  in  Chapter  VI.  It  includes  comparisons  of 
spectra;  radial,  axial,  and  azimuthal  group— flux  profiles; 
and  tritium  production  rates  in  the  different  blanket 
models.  Advantages  of  constructing  the  FRBF  at  Purdue  and 
of  the  use  of  this  facility  to  produce  more  and  better 
fusion  reactor  blanket  neutronics  and  tritium  breeding 
experimental  data  are  discussed  in  chapter  VII.  This 
chapter  discusses  safety  issues,  predicted  measurement 
precision  and  accuracy,  applicability  of  this  design  to  the 
more  complicated  geometries  of  fusion  reactors,  and 
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flexibility  for  multiple  and  inexpensive  benchmark 
experiments.  Finally,  the  research  is  summarized  in  Chapter 
VIII,  and  recommendations  for  further  research  are 
presented. 


II.  REVIEW  OF  FUSION  REACTOR  BLANKET  EXPERIMENTS 


Few  actual  fusion  blanket  experiments  have  been 
conducted.  Leonard  reviewed  the  blanket  experiments  up  to 
19752 ,  and  Maynard  reviewed  essentially  the  same  experiments 
at  the  BNL  cross  sections  meeting  in  1977. 3  L.  Green  then 
reviewed  the  experiments  reported  up  to  1980. 4  In  1984 
Woodruff  summarized  Green's  paper  and  reviewed  the  fusion 
blanket  experiments  up  to  that  time.5  These  experiments  are 
summarized  in  Table  II.  1,  and  the  corresponding 
computational  methods  are  summarized  in  Table  II. 2.  All  the 
experiments  except  LOTUS  and  the  LBM  used  a  neutron  source 
produced  by  directing  a  deuteron  beam  on  a  metal  (Cu  or  Ti) 
target  with  tritium  imbedded  in  it.  However,  different 
researchers  treat  the  source  differently  in  their 
measurements  and  calculations.  Measurements  have  included 
local  tritium  production  rates  (Tn  is  the  tritium  production 
rate  in  natural  lithium,  tr itons/cra3-source  neutron,  T6  from 
®Li ,  T7  from  7Li),  neutron  and  gamma-ray  spectra  (including 
angular  flux),  and  fission  rates.  A  brief  discussion  of  the 
reported  experiments,  in  chronological  order,  follows  the 


Table  II.  1 


A  Summary  of  Tritium  Breeding  Experiments 


Ref. 

Year 

Place 

Geometry 

Materials 

6 

1954 

LANL 

60  cm  sphere 

nLiD 

7-9 

1974-77 

Karlsruhe 

1  m  sphere 

nLi ,  SS 

10&11 

1974-75 

Julich 

120  cm  cyl 

nLi ,  SS 

12 

1976 

Jul ich 

120  cm  cyl 

"Li.  SS 

13-17 

1975-78 

JAERI 

pseudo— sphere 

"Li,  U,  C. 
"Li  a0 

18 

1978 

JAERI 

pseudo— sphere 

"Li aO ,  C 

19 

1978 

LANL 

60  cm  sphere 

•LiD 

20 

1979 

Osaka  U. 

slabs 

"Li  or  C 

21 

1980 

Tokyo  U. 

slab 

"LiF 

22 

1981 

U.  Birin.  , 
England 

1.25  m  sphere 

"LiF  powder 
A1 

23 

1981 

ORNL 

slabs 

LiH,  Pb-LiH 

24 

1980-81 

UCSB 

88  cm  cyl 

"Li,  SS 

25 

1983 

JAERI 

pseudo-cyl  slab 

"LiaO 

26 

1983 

JAERI 

pseudo— cyl  slab 

"Li  a0 

27-29 

1984—? 

Lausanne 

LOTUS 

Fusion- 

fission 

30 


1986-? 


PPPL-TFTR 


LBM 


"Lia0,  SS 
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Table  II. 2.  Computational  Methods  in  Tritium  Breeding 
Experiments . 


Ref. 

Codes 

Cross 

Sections 

Source 

Treatment 

6 

DTF-IV  (Sn) 

ENDF/B— I I I 

isotropic 

7-9 

DTK  (Sn) 

ENDF/B— I I I 
KEDAK3  (Fe) 

meas'd  spectrum 
anisotropic 

10&11 

ANISN,  DOT— I I 

MORSE 

ENDF/B— I I I 

same  as  7 

12 

ANISN,  MORSE 

ENDF/B— I I I 

isotropic,  1  gp 

DOT- I I 

13-17 

ANISN 

ENDF/B— I I I 

isotropic,  1  gp 

18 

ANISN 

ENDF/B- IV 

isotropic 

19 

MCN 

(Monte  Carlo) 

ENDF/B— III 

67  UK-LASL 

E/0  correlation 
anisotropic 

20 

ANISN.  NITRAN 

NIMOS 

ENDF/B- IV 

measured , 

14.8  MeV 

21 

ANISN 

ENDF/B— I I I&IV 

measured , 

14  MeV 

22 

MORSE.  ANISN 

ENDF/B- IV 

anisotropic , 
meas'd  spectrum 

23 

DOT 

ENDF/B- IV 

E/0  correlation 
anisotropic 

24 

MORSE-L 

ENDF/B- IV 

disc,  1  group, 
isotropic 

25 

BERMUDA-2DN 

ENDF/B- IV 

meas'd  spectrum 

26 

same  as  25 

27-29 

ANISN,  MCNP 

DLC-2 ,  SINEX 

anisotropic , 
col lided  MCNP 

30 

Monte  Carlo 

p 1 asma 
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The  first  documented  tritium  breeding  integral  experiment 
was  conducted  at  LANL  in  1954,  but  the  neutronics  analysis 
wasn't  performed  until  1972  by  Muir  and  Wyman.*  The  60-cm— 
diameter  natural  lithium  deuteride  (nLiD)  sphere  with  a 
neutron  source  at  the  center  produced  T7  and  Tn  data.  These 
data  were  compared  to  1-dimensional  breeding  predictions, 
and  discrepancies  ranged  from  a  few  to  ~35X,  depending  upon 
isotope  and  radius.  They  argued  that  the  7Li(n,n'a)T  cross 
sections  were  too  great,  but  not  enough  to  explain  the 
magnitudes  of  errors  in  the  data. 

Another  experiment  was  conducted  in  Great  Britain  in 
1964,  but  information  on  this  experiment  is  difficult  to 
obtain  (report  AWRE  NR-4/64  is  supposed  to  have  information 
on  this  experiment).  Maynard  listed  this  cylindrical  nLi 
experiment  in  his  table  in  Ref.  3,  but  gave  no  comparison  of 
results  to  calculations. 

Various  transport  codes  were  compared  to  the  results  of 
experiments  in  a  1-m-diameter  nLi  sphere  with  a  stainless 
steel  shell  from  1974  through  1977. 7"*  Bachmann  et  al . 
found  that  Tn  was  overpredicted  by  ~35%  throughout  the 
blanket,  and  neutron  flux  was  greatly  overpredicted  between 
3  and  11  MeV.  The  neutron  source  in  this  experiment  was 
strongly  anisotropic,  so  it  was  measured  and  the  measured 
spectrum  was  used  as  the  source  in  their  codes. 


Herzing  et  al .  performed  experiments  in  a  cylindrical 
blanket  model  of  nLi  ab  Julich  in  1975. 10,11  They  measured 
T7,  T# ,  and  Tn  and  speebra  using  solid  sbabe  brack 
debecbors.  They  used  various  1.  2,  and  3-dimensional  codes 
and  assumed  an  isobropic  neubron  source.  They  found  bhab 
bhe  resulbs  of  Monbe  Carlo  predicbions  mabched  bhe  Tn  daba 
excepb  ab  bhe  boundaries.  Since  obhers  have  debermined  bhab 
bhe . 7Li (n ,n 'a)T  cross  secbions  are  ~15X  boo  large  and  bhe 
scabbering  cross  secbions  are  quibe  inaccurabe  in  ENDF/B- 
III . s  ,  •  ,  7  they  should  have  found  discrepancies. 

Addibional 1 y ,  bhey  debermined  T7  was  overpredicbed 
bhroughoub  bhe  blankeb  by  aboub  15-20%.  Only  bhe  Monbe 
Carlo  mebhod  gave  bhem  whab  bhey  considered  accepbable 
resulbs.  This  experimenb  was  repeabed  in  1970  wibh  a  15-cm- 
bhick  graphibe  blankeb  added.10  Again  T* ,  T, ,  and  T7  were 
measured.  Generally  good  agreemenb  was  obbained  in  bhe 
ouber  region  of  bhe  blankeb,  bub  Tn  was  underpredicbed 
toward  bhe  cenber  (~12%)  and  ab  bhe  ouber  regions  of  bhe 
reflecbor  (~16%) ,  and  overpredicbed  wibhin  5  cm  of  bhe 
blankeb/ref lecbor  inberface  (~8%) .  Additionally,  T,  was 
greably  overpredicted  (~45%)  ab  the  inner  wall.  Again, 
ENDF/B-III  cross  secbions  were  used,  and  T7  should  have  been 
overpredicbed. 

A  series  of  experiments  were  conducted  ab  JAERI  between 
1975  and  1979.ia-1#  Pseudo— spherica 1  blankets  were  formed 
by  assembling  rectangular  blocks  of  various  materials  (Li, 


U,  C,  LiaO).  Fission  rates  were  measured  in  miniature 
fission  chambers  containing  aa#U,  aa#U,  aaaTh,  and  aa7Np. 
Discrepancies  between  ANISN  calculations  and  measurements 
depended  on  blanket  material,  fission  isotope,  and  depth; 
they  ranged  from  +15*  to  -15*.  sometimes  with  strong 
transitions  through  the  blanket.  Angular  flux  spectra  were 
compared  to  predictions  in  the  last  experiment;  results 
indicated  that  C/E  was  a  strong  function  of  angle,  with 
differences  up  to  30*. 17  Maekawa  et  al.  attributed 
discrepancies  to  inadequate  cross  section  libraries  and  to 
non-spherical  interfaces  due  to  the  use  of  the  blocks. 

Another  60— cm-diameter  spherical  assembly  was  tested  at 
LANL  in  1978. 19  It  contained  »LiD  (96.6*  *Li).  and  breeding 
was  measured.  Tc  and  T7  were  measured  and  compared  to  3- 
dimensional  Monte  Carlo  predictions  (using  ENDF/B-III).  The 
source  was  treated  with  an  angle— energy  correlation  and 
anisotropic  yield.  T»  was  generally  very  close  to  the 

••  f 

predictions  (except  at  boundaries),  but  T7  was  again 
overpredicted  by  about  10*  (except  at  the  outer  boundary). 

Yamamoto  et  al.  measured  angular  spectra  from  various 
thicknesses  of  Li  or  graphite  slabs  in  1979. 20  They 
compared  the  measurements  with  various  codes,  using  a 
monochromatic  14.8  MeV  source.  Transverse  leakage  was 
blamed  for  a  large  part  of  the  20-30*  discrepancies  in  the 
spectra  at  low  energies.  The  spectra  were  much  softer  than 
predicted,  especially  below  8  MeV  for  the  Li  slabs.  ANISN 
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results  were  totally  unacceptable. 

Another  slab,  this  time  of  natural  LiF,  was  tested  at  the 
University  of  Tokyo  in  1980. 21  Neutron  spectra  and  Tn  were 
measured  and  compared  to  ANISN  predictions,  which  were  found 
to  be  good  to  no  better  than  20-302. 

LiF  was  also  tested  in  a  spherical  blanket  model  in 
1980. 22  However,  in  this  experiment  a  1 . 25— m— diameter 
aluminum  sphere  was  filled  with  a  powder — variable  packing 
density  may  have  affected  the  results.  An  anisotropic, 
measured  spectrum  was  used  as  the  source  in  various  codes. 
Great  attention  to  detail  in  this  experiment  did  not  seem  to 
ensure  better  results.  Some  measured  spectra  were  much 
softer  than  predicted,  others  were  much  harder.  Woodruff 
was  concerned  with  this  because  of  the  quality  of  the  effort 
and  results  which  he  considered  opposite  to  those  of  Profio 
et  al .  (see  Ref.  24  below). 

Experiments  were  conducted  on  slabs  of  LiH  with  and 
without  5  cm  of  Pb  at  ORNL  in  l'981.2a  Santoro  et  al. 
measured  neutron  and  Y-ray  leakage  spectra  and  compared 
their  results  to  DOT  predictions.  They  used  a  first 
collision  source  produced  by  a  code  which  used  an  angle— 
energy-yield  correlation  for  neutron  production  in  the 
target.  Results  were  good  between  6  and  10  MeV,  but  spectra 
were  underpredicted  above  and  below  that  range.  This  may  be 
due  to  an  error  in  the  source  correlation  in  the  forward 


direction  (toward  the  slab).83 


Profio  et  al.  experimented  with  a  cylindrical  nLi  blanket 
at  UCSB  in  1980. 34  This  system  was  unique  due  to  the  radial 
rather  than  axial  orientation  of  the  d— beam  tube.  They 
measured  fast  scalar  flux  alone  the  axis  of  the  cylinder, 
rather  than  radially.  Monte  Carlo  calculations  using  ENDF/ 
B— IV  significantly  underpredicted  the  penetration  of  fast 
neutrons,  by  50-100*  for  some  spectral  peaks. 

Maekawa  et  al .  measured  the  leakage  spectra  from  pseudo— 
cylindrical  slabs  of  Lia0  by  t ime-of-f 1 ight  methods  in 
1983. 38,38  They  used  a  2-dimensional  transport  code  with 
ENDF/B-IV  cross  sections  to  predict  the  high-energy  angular 
neutron  flux.  Measurement  errors  were  1-20%  depending  on 
slab  thickness  and  angle.  Spectral  shapes  were  reproduced, 
but  calculated  spectra  were  softer  than  measured  spectra, 
and  some  spectral  peaks  were  overpredicted  by  ~50%.  Aside 
from  the  past  experiments  which  were  previously  discussed, 
and  despite  the  obvious  requirement  for  further  research, 
only  two  new  projects  are  currently  in  progress. 

One  new  project  for  testing  conceptual  designs  for  fusion 
and  fusion-fission  hybrid  blankets  is  being  conducted  at 
Lausanne,  Switzer land. 37  It  is  called  the  LOTUS  project  and 
uses  a  higher  intensity  neutron  source  than  most  other 
blanket  experiments.  Sitaraman  et  al.  did  a  Monte  Carlo 
analysis  of  its  HAEFELY  mixed-beam  DT— neutron  generator.38 
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They  determined  that  the  neutron  source  for  LOTUS  would  be 
highly  anisotropic,  and  noted  only  37%  of  the  neutrons  would 
be  emitted  toward  the  blanket;  of  those,  only  67%  would  be 
uncollided  neutrons.  Initial  experiments  were  scheduled  for 
1985,  and  the  neutronics  predictions  have  been  completed. 
Sahin  and  Kumar  have  predicted  that  lateral  leakage  will  be 
a  major  source  of  discrepancies."  Experimental  results 
have  not  appeared  in  the  open  literature. 

Another  planned  experiment  will  be  the  Lithium  Blanket 
Module  (LBM)  to  be  placed  in  the  TFTR  at  PPPL.  Jassby  et 
al.  describe  the  program,  including  design,  manufacturing 
techniques,  neutronics  analyses,  and  experiment  plans.*0 
The  LBM  will  be  a  segment  81  cm  x  76.7  cm  x  86.5  cm,  with 
LiaO  pellets  contained  in  steel  rods  and  with  a  steel 
reflector.  Tritium  production  and  neutron  current  are 
predicted  to  be  flat  across  the  LBM  to  within  ~5%  except 
within  8  cm  of  the  sides  (probably  due  to  neutron  leakage 
from  the  surrounding  structure).  When  the  TFTR  burns 
deuterium  fuel  within  the  next  few  years,  the  fusion 
community  will  have  its  first  real  fusion  reactor  blanket 
data. 

The  above  experiments  suffer  from  many  deficiencies: 
complicated  geometries,  problems  with  anisotropy  of  the  DT— 
neutron  source,  inadequate  measurement  techniques,  and 
designs  which  were  nob  flexible.  Only  the  JAERI ,  ORNL,  and 
LOTUS  experiments  allow  changes  in  blanket  composition;  the 


ORNL  system  and  LOTUS  are  thick  slabs,  and  JAERI  ussd  non- 
smooth  boundaries.  Thus,  many  of  Leonard's,  Maynard’s, 
Green’s,  and  Woodruff’s  conclusions  are  still  valid:  T7 
cannot  be  predicted  accurately,  measurement  techniques  need 
refinement  in  both  accuracy  and  sensitivity,  cross  section 
libraries  and  transport  codes  need  improvements,  and 
integral  experiments  are  needed  on  prototypic  fusion  reactor 
blanket  designs  to  validate  the  codes  and  the  designs.  The 
lack  of  fusion  reactor  blanket  experiments  since  1984  should 
cause  concern,  since  the  past  research  has  not  yielded  data 
which  validated  transport  codes  or  cross  sections  for 
tritium  breeding  and  neutronics  predictions. 

The  deficiencies  previously  discussed  have  contributed  to 
the  deviations  between  experimental  results  and  theory 
sufficiently  to  prevent  validation  (or  invalidation)  of 
neutronics  codes.  A  new  fusion  reactor  blanket  simulation 
facility  without  the  deficiencies  in  design  of  the  past 
experiments  will  provide  data  which  can  be  used  to 
inexpensively  validate  those  codes  and  the  cross  section 
libraries  for  a  variety  of  fusion  reactor  blankets.  The 
following  chapter  describes  how  the  FBBF  was  designed  for 
this  same  purpose  (for  fast  reactor  neutronics),  and  it 
summarizes  the  highly  successful  results. 


III.  OVERVIEW  OF  THE  FBBF 


III. 1 .  Facility  Description 

Ott,  Clikeman.  and  Harms  describe  the  FBBF  in  Ref.  31, 
and  present  a  summary  and  analysis  of  the  many  experiments 
conducted  in  it.  The  facility  includes  a  heavily  shielded 
room,  which  contains  the  cylindrical  blanket  simulation 
system,  and  an  exterior  control  and  instrumentation  room. 

The  experimental  system  itself  consists  of  an  axial  neutron 
source  of  assCf  of  about  109  n/s ,  a  two-region  spectrum 
modifier  with  an  outer  radius  of  22  cm,  a  breeding  blanket 
with  an  outer  radius  of  73  cm  (51  cm  thick),  and  a  NaCl  and 
SS>  reflector  with  an  outer  radius  of  89  cm  (16  cm  thick). 
Four  pellets  of  Cf  are  distributed  on  the  axis  to  create  the 
desired  neutron  source;  they  are  fixed  in  a  moveable  rod 
(lowered  into  underground  shielding  while  changing 
experiments,  etc.),  and  their  intensities  and  locations  are 
such  that  a  cosine-shaped  source  results  after  some 
distance.  The  inner  region  contains  4 . 835-enr  iched ,  close 
packed  U0a  fuel  rods  clad  in  SS ,  with  the  gaps  between  rods 
filled  with  B«C  powder;  the  outer  modifier  contains  the  same 
fuel  rods,  but  with  wider  spacing,  and  the  gaps  contain  Na- 
filled  SS  tubes.  These  regions  modify  the  neutron  spectrum 
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from  the  Cf  source  to  approximate  the  leakage  spectrum 
entering  the  blanket  from  a  typical  1000— MWe  fast  reactor 
core.  They  also  smooth  the  flux  from  the  four  sources  so 
that  the  axial  flux  distribution  at  the  modifier/blanket 
interface  closely  matches  that  which  would  be  produced  by  a 
continuously  distributed  chopped-cos ine-shaped  source  (as  in 
a  critical  fast  reactor). 

The  blanket  region  contains  natural  U02  fuel  pellets  clad 
with  Al,  then  placed  in  A1  or  SS  secondary  cladding.  The 
rods  are  placed  in  a  uniform  hexagonal  array  with  the 
correct  pitch  between  fuel  rods,  and  no  other  materials  are 
used  (voids  contain  only  air).  All  experiments  are  placed 
between  fuel  pellets  or  in  a  vacancy  created  by  removing  an 
experimental  rod.  Thus,  measurements  can  be  made  throughout 
the  blanket  in  axial,  azimuthal,  or  radial  traverses  with 
negligible  perturbations  of  the  flux.  Since  Al  has  about 
the  same  macroscopic  cross  section  and  slowing  down  power  as 
Na,  it  plays  the  role  of  the  Na  coolant  in  a  fast  reactor. 
These  features  allow  the  FBBF  to  simulate  very  closely  the 
blanket  of  a  full-scale  fast  breeder  reactor — in  geometry, 
in  source  spectrum  and  current,  and  in  materials. 

A  new  blanket  is  now  being  constructed  with  the  fuel  rods 
arranged  in  hexagonal  arrays  (assemblies).  Additionally, 
initial  plans  are  being  formulated  for  fabricating  several 
thousand  metal  fuel  rods  for  use  in  the  FBBF  for  testing  the 
neutronic  characteristics  of  the  Integral  Fast  Reactor.  The 


outstanding  results  obtained  in  past  F3BF  experiments, 
discussed  in  the  following  section,  should  also  be  found 
with  these  new  blanket  configurations. 

III. 2.  Results  of  FBBF  Experiments. 

Measurements  in  the  FBBF  have  included  neutron  capture 
rates  in  2  3  8tJ  and  232Th;  fission  rates  in  235U,  23  8U,  232Th, 
and  239Pu  by  solid  state  track  detector  (SSTD)  measurements; 
resonance  neutron  capture  rates  in  286W,  S5Mn,  and  197Au; 
helium—  and  croton— reco i 1  neutron  spectrometry;  and  gamma- 
ray  energy  deposition  rates  in  thermo- 1 uminescent  dosimeters 
(TLD).  These  measurements  have  been  compared  to  diffusion 
calculations  based  on  standard  50-group  cross  sections.  The 
results  of  these  comparisons  are  expressed  as  ratios  of 
calculation  to  experiment  (C/E  values,  C  is  based  on 
absolute  neutron-source  strength,  E  on  absolute 
measurements).  The  general  trend  in  these  C/E  values  is  to 
begin  at  ~1  at  the  transformer/blanket  interface,  then  to 
decrease  linearly  with  increasing  radius  to  ~0 . 65  at  the 
outer  circumference  of  the  blanket.  The  cause  of  this  trend 
has  not  been  determined,  but  is  suspected  to  be  due  to  a 
gross  (bulk)  neutron  effect,  which  should  appear  in  all  C/E 
results.  Exeptions  to  this  trend  should  be  examined 
carefully  for  errors  in  measurement  technique  or 
interpretation,  or  in  neutronics  calculation.  This  careful 
examination  has  led  to  the  identification  of  many  problems 
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in  blanket  neutronios:  faulty  measurement  techniques,  errors 
in  treating  resonance  self-shielding  of  detection  isotopes 
in  group  constant  production  codes,  shielding  of  resonance 
absorption  in  low-concentration  materials  by  high- 
concentration  materials,  and  short-range  deformations  of 
reaction  rates  (due  to  spectrum  transition  and  its  effect  on 
se 1 f-shie 1  ding) . 

The  consistency  of  the  results  of  the  various  measurement 
techniques  employed  in  the  FBBF  experiments  lends  great 
credibility  to  the  results.  The  careful  design  of  the  FBBF 
insured  these  outstanding  results.  The  conclusion  of  Ott  et 
al.  in  Ref.  31  that  some  of  the  trends  should  be  relevant  to 
fusion  reactor  blankets  came  from  a  careful  analysis  and 
synthesis  of  the  varied  measurement  results.  They  predicted 
that  one  result  will  be  especially  applicable  to  fusion 
reactor  neutron  transport:  the  difficulties  in  describing 
the  transitory  neutron  fluxes  in  the  FBBF  blanket  should  be 
more  severe,  due  to  the  wider  range  of  neutron  energy  in  a 
fusion  blanket.  Thus,  the  idea  for  the  FRBF  was  born 
naturally  from  the  results  of  FBBF  investigations.  However, 
the  feasibility  of  constructing  the  FRBF  and  using  it  to 
produce  more  and  better  test  data  remains  to  be 
demonstrated:  that  demonstration  follows  in  the  next  four 
chapters  of  this  thesis. 


IV.  CONCEPTUAL  DESIGN  OF  A  FUSION  REACTOR  BLANKET  FACILITY 


IV. 1.  Objectives 

The  general  objectives  of  the  research  reported  herein 
were  threefold: 

l.a)  to  demonstrate  the  advantages  of  constructing  a 
fusion  reactor  blanket  testing  facility  (the  FRBF)  which  is 
an  analog  of  the  FBBF ,  and 

b)  to  show  the  feasibility  of  constructing  the  FRBF 
at  Purdue  University: 

2)  to  demonstrate  the  new  facility  will  produce  more 
and  better  data  and  results  than  are  (or  will  be)  available 
elsewhere;  and 

3)  to  show  the  results  are  applicable  to  fusion 
reactors  although  neither  the  geometry  nor  the  source  are 
exactly  duplicated. 
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IV. 2.  General  Considerations  for  the  Design 


Before  beginning  the  design  of  the  FRBF,  the  requirements 


of  a  fusion  reactor  blanket  simulation  facility  must  be 


analyzed.  Then,  convincing  arguments  must  be  made  which 


prove  that  the  FRBF  can  be  built  as  an  FBBF-analog  which 


will  satisfy  these  requirements.  In  Chapter  II  the  multiple 
problems  with  past  experiments,  current  facilities,  and 
planned  experiments  were  demonstrated.  Those  problems  can 
be  summarized  as: 

•  severe  source  anisotropy, 

•  inappropriate  geometry, 

•  ray  effects  due  to  point  sources, 

•  short  measurement  times  due  to  finite  lifetime  of 
targets,  and 

•  inflexible  designs. 

To  show  that  experiments  conducted  in  the  FRBF  will 
contribute  to  the  advancement  of  fusion  reactor  technology, 
this  research  needs  only  to  demonstrate  that  and  how  those 
problems  are  eliminated,  reduced,  or  ameliorated.  Although 
the  blanket  model  is  neither  toroidal ly  shaped  (as  in  a 
Tokamak)  nor  semi— inf ini te 1 y  long  (as  in  a  large  mirror), 
nor  is  the  source  a  "tube"  of  neutron— producing  reactions, 
the  results  must  be  applicable  to  those  geometries.  Once 
precise  data  are  produced  which  lead  to  improved  transport 
codes  and  validation  of  data  libraries,  then  confidence  in 
these  codes  to  predict  transport  in  real  blankets  will  have 
been  obtained. 

The  results  of  the  research  reported  in  this  paper  will 
simultaneously  satisfy  the  three  objectives  previously 
mentioned.  The  remainder  of  this  chapter  will  present  the 
desired  properties  of  a  fusion  reactor  blanket  experiment. 


and  then  will  describe  how  these  properties  can  be  supplied 
by  the  appropriate  design. 

IV. 3.  Desired  Features  of  a  FRB  Simulation  Facility 

The  facility  for  fusion  reactor  blanket  simulation 
experiments  must  have  certain  properties  so  that  it  will 
produce  improved  data.  This  new  data  must  be  more 
meaningful  than  that  which  has  been  produced  in  the  past;  or 


it  must 

be  data 

which 

might 

not  otherwise 

be  produced  in  the 

future , 

due  to 

the  hi 

gh  cost 

of  one— of— a— 

kind  simulations  or 

fusion 

reactor 

exper i 

ments . 

The  required 

properties 

inc lude ; 

•  simple  geometry  and  an  easily  modelled  source  for 
improved  analysis; 

•  flexibility  for  changing  blanket  materials  or 
composition  at  moderate  cost; 

•  long  measurement  times  for  greater  accuracy  and  • 
precision ; 

•  multiple  measurement  techniques  for  verification  of 
results  and  broader  comparisons; 

•  low  cost  of  changes  in  composition  or  blanket 
structure  for  many  experiments;  and 

•  low  cost  of  initial  construction  and  operation. 
These  properties  are  synergistic — geometry  affects  analysis, 
cost,  and  flexibility;  source  affects  analysis  and 
measurement  statistics;  cost  affects  number  and  variety  of 
experiments  and  analyses;  etc.  Therefore,  proper  selection 


of  the  geometry  and  source  can  lead  to  an  experimental 


system  which  can  yield  better  data  quickly,  accurately,  and 
cheaply. 

IV. 4.  Scope  of  Research 

The  research  which  needed  completion  prior  to  a  decision 
to  proceed  with  engineering  design  of  the  FRBF  and  its 
source  was  to  develop  a  preliminary  conceptual  design;  then 
to  analyze  that  design  to  determine  if  it  meets  the 
guideline  previously  set  forth.  That  is,  does  it  have  a 
geometry  and  a  source  that  can  be  analyzed  by  standard 
methods  without  significant  deviations  due  to  inadequate 
representabi 1 ihy?  The  steps  to  be  followed  were  set  forth 
as : 

1)  select  a  blanket  design, 

2)  select  a  source  arrangement,  and 

3)  perform  neutronics  analyses. 

These  three  efforts  are  discussed  thoroughly  in  the 
following  three  sections. 

IV. 5.  Blanket  Design 

The  overall  blanket  geometry  for  the  FRBF  has  already 
been  selected  due  to  the  success  of  the  FBBF.  The  non¬ 
availability  of  appropriate  neutron  sources  for  slab  (plane 
wave)  or  spherical  (isotropic)  geometries  also  leads  to  the 
selection  of  cylindrical  geometry,  as  does  the  requirement 
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for  flexibility.  Additionally,  others  have  performed 
optimization  studies  which  have  shown  the  FBBF  blanket  is 
about  the  correct  size  for  fusion  blanket  research.10"12 
The  final  blanket  design  includes  an  inner  source  region  and 
vacuum  tube,  a  first  wall  of  variable  thickness,  a  blanket 
region,  a  reflector  or  plenum  region,  and  a  steel  shield. 

The  blanket  designs  are  based  on  both  the  design  of  the 
present  FBBF  and  on  the  primary  candidates  selected  in  the 
Blanket  Comparison  and  Selection  Study  (BCCS) ,  which  was 
managed  by  ANL  for  the  U.S.  Fusion  Power  Program  (FP P) . 1 
The  thicknesses  of  blanket  regions  are  the  order  of  those 
reported  for  fusion  reactor  models  (based  on  Tokamak  or 
large  mirror  reactors)  in  the  FINESSE  study:  but  the  radius 
of  the  plasma  region  is  reduced  by  about  an  order  of 
magnitude,  and  the  reflector/shield  region  is  decreased  in 
thickness.  The  FRBF  will  actually  be  about  the  same 
dimensions  as  the  proposed  fusion  test  system,  also 
described  in  FINESSE.  This  system  will  be  placed  in  the 
center  of  the  MFTF—B  at  LLNL,  where  tremendous  magnetic 
fields  will  confine  the  test  plasma  inside  the  test  blanket 
(at  a  probable  cost  of  millions  of  dollars  per  test).  Of 
the  blanket  candidates  selected  by  the  BOSS,  only  two  were 
selected  for  analysis  in  this  study. 

Four  blankets  were  selected  as  primary  candidates  by  the 
BOSS  group;  two  of  those  were  selected  for  this  research. 
They  are: 


•  Blankeb  HeC — a  helium-cooled  solid  blanket  design, 
with  natural  lithium  composition;  and 

•  Blanket  LLC — a  1 iquid-meta 1 -coo  1 ed  liquid  blanket 
design,  with  neutron  multiplication  and  enriched  lithium. 

These  two  models  offer  totally  different  neutronic 

environments  for  studying  the  effects  of  asymmetries  and 

anisotropic  sources;  thus  helping  insure  that  the 

conclusions  reached  from  comparing  the  results  of 

computations  are  not  blanket-design  or  composition 

dependent.  The  general  FRBF  blanket  geometric  design  is 

shown  in  Figures  IV. 1. A  and  IV.l.B.  One  dimensional 

representations  of  the  complete  blanket  designs,  including 

materials  and  densities,  are  shown  in  Figures  IV. 2. A  and 

IV. 2. B.  The  blanket  inner  radius  and  outer  radius  are  22 

and  100  cm  respectively,  the  height  of  the  blanket  is  100 

cm,  and  it  is  capped  by  a  20  cm  thick  steel  reflector/ 

shield.  Number  densities  of  the  components  of  the  blanket 

materials  are  listed  in  Appendix  A.  The  blanket  will  be 

constructed  as  a  hollow  annulus,  with  removable  azimuthal 

and  axial  segments  to  allow  flexibility  in  design  of 

experiments.  The  source  which  will  be  created  in  the  center 

of  the  blanket  is  a  critical  element;  its  conceptual  design 

is  described  in  the  next  section. 

IV. 6.  Neutron  Sources 

A  neutron  source  arrangement  for  the  FRBF  was 
investigated  by  Aparcedo  in  an  earlier  paper.33  He  studied 
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Figure  IV. 2. A.  1-Dimensional  Model  of  Blanket  HeC. 
cooled  blanket  with  steel  structure. 
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Figure  IV. 2. B.  1-Dimensional  Model  of  Blanket  LLC. 

Lii7Pb33  self-cooled  with  Primary  Candidate  Alloy  structure. 


a  multi-target  arrangement  along  the  axis  of  a  cylinder 
which  would  produce  an  approximately  cosine— shaped  total- 
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flux  distribution  at  the  inner  wall.  In  this  preliminary 
investigation,  he  used  an  isotropic  emission  of  T(d,n)a 
neutrons  from  flat  annular  TiT  targets;  the  result  was  a 
distribution  of  point  sources  which  was  symmetric  about  the 
mid— plane  of  the  cylinder.  In  fact,  the  neutrons  emitted 
from  such  targets  are  anisotropic  in  both  energy  and 
yield. 32 

A  kinematic  treatment  of  the  two-body  reaction,  using 
conservation  of  momentum  and  energy  in  the  centei — of— mass 
(COM)  reference  frame,  yields  the  energy-angle  relationship 
for  the  T(d,n)a  reaction: 

En(fl)  -  BET (cosS  +  (D/B  -  sin2©)1'2)2,  (4.1) 

where 

En(e)  »  energy  of  the  neutron  at  the  angle  Q 
from  the  d— beam, 

Et  «  Ed  +  Q  (Q  *  17.59  MeV),  and 

Ed  ■  d-beam  energy. 

The  solid  angle  ratio  for  conversion  of  the  COM  reaction 
cross  section  to  the  laboratory  cross  section  is 

Sf(S)  -  En(e)/ET(AC)l/a(D/B— sina0)1/a.  (4.2) 

The  dimensionless  constants  A,  B,  C,  and  D  are  determined 
from  the  masses  and  energies  of  the  interacting  species. 

For  the  T(d,n)a  reaction,  they  are  given  by 
A  -  MrMdMa(Ed/ET). 

B  -  MrMdMn(Ed/ET) . 


where 


Hr  «  l/(Md  +  Mt)(Mn  +  Ma).  Note  that  the 

product  of  nr  with  two  masses  is  dimensionless  also. 

When  equations  4.1  and  4.2  are  used  to  determine  the 
energy  and  yield  of  neutrons  emitted  from  a  deuteron 
incident  on  a  TiT  target,  the  result  is  an  angle-energy 
correlation  of  neutron  emission.  Table  IV. 1  lists  the 
angle— energy  correlation  for  a  150  KeV  deuteron. 


Table  IV. 1.  Angular  Distribution  of  d-T  Neutrons. 


The  energy  of  neutrons  and  alpha  particles  in  the  table  is 
in  MeV.  Note  that  the  sum  of  En  and  Ea  is  17.74  MeV,  which 
is  the  sum  of  the  reaction  energy,  Q,  and  the  deuteron 
energy,  Ed ,  at  each  angle.  Based  on  the  solid  angle  ratios 
listed,  the  neutron  production  rate  is  about  12%  greater  in 


the  forward  direction  than  in  the  backward  direction. 

The  treatment  which  determines  this  distribution  is 
incomplete  for  a  deuteron  beam,  however.  Due  to  slowing 
down  and  scattering  of  the  beam  in  the  target  material ,  the 
reaction  energy  is  not  equal  to  the  beam  energy  for  most 
deuterons;  in  fact,  the  majority  of  reactions  occur  at  the 
resonance  energy  of  109  KeV  for  the  T(d,n)a  reaction. 
Therefore,  the  angle— energy  correlation  must  be  determined 
as  a  function  of  the  slowing  down;  it  has  been  evaluated  and 
reported  with  respect  to  this  loss  of  energy  before 
reaction. 32  The  implementation  of  this  angl e— energy-s 1  owing 
down  relationship  will  be  presented  in  Chapter  V.  Since 
there  is  a  distribution  of  energy  and  yield  from  the  d— t— 
neutron  source,  and  because  of  concerns  of  tritium  leakage 
from  the  system,  other  high-energy  neutron  sources  were 
considered. 

Other  reactions  which  yield  high-energy  neutrons  were 
also  evaluated,  e.g.,  the  d— 6Li  and  d— 7Be  reactions  were 
considered.  These  reactions  have  even  stronger  variations 
in  energy  and  yield  than  the  d-T  reaction,  because  they 
require  extremely  high  beam  energy  for  acceptable  neutron 
production  rates.  For  example,  the  energy  and  solid  angle 
ratios  for  the  eLi(d,n)7Be  reaction  (Q  ■  3.38  MeV)  for  a  15 
MeV  deuteron  are  shown  below  in  Table  IV. 2.  Neutron  energy 
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2iZi<t>iSz£  -  0.  (4.3) 

where  z{  is  the  axial  position  at  the  first  wall  (:  ■  0  at 
the  mid-plane),  <$>£  is  the  flux  (a  function  of  Zj  and  Sj  and 
of  z£),  and  all  Sz£'s  are  equal.  The  result  of  this  search 
for  the  central  source  position  was  Z3  *  —1.2  cm  (offset 
toward  the  beam).  Two  additional  source  positions  (Z2  and 
Z4)  were  added,  and  another  constraint  was  imposed: 

,  j  5zi  =  °»  and  (4.4) 

AZ2  =  AZ3,  (4.5) 

where  AZj  =®  Z j «.  £  —  Zj  ,  the  separation  between  the  sources. 
The  result  of  this  search  showed  that  the  position  of  the 
central  source,  Z3 ,  was  very  weakly  dependent  upon  the 
source  separation,  AZ  (±0.02%  over  a  range  of  AZ  *  0  to  30 
cm)  . 

To  find  the  optimum  positions  for  three  sources,  a  least- 
squares-error  routine  was  used  to  minimize  the  difference 
between  the  flux  profile  and  the  chopped— cosine  function. 

The  error  routine  minimized  e2,  where 

e3  -  2i  (4»i-cos(z£*)  )2cos2(z£*)  ,  and  (4.6) 

z  £  *  »  rrz  £  /H  ’  . 

H',  the  extrapolated  height,  was  based  upon  the  reflector 
savings  of  the  FBBF.  The  summation  interval  ignored  10  cm 
at  each  end  of  the  cylinder,  since  slight  errors  there  will 
have  a  minimal  effect  at  the  midplane,  where  measurements 
will  be  concentrated.  The  extra  cos2  term  in  equation  4.6 
is  a  weighting  factor  which  tends  to  improve  the  fit  around 
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the  mid— plane.  The  results  of  this  search  are  illustrated 


in  Figure  IV. 3  (along  with  the  flux  profile  produced  by  a 


single  anisotropic  point  source).  The  y-coordinate  of  the 
octagons  represents  the  relative  strengths  of  the  three 
sources.  The  RMS  error  over  the  central  80%  of  the  cylinder 
is  about  1%.  Although  it  appears  that  three  sources  may  be 
adequate,  considering  the  projected  accuracy  of 
measurements,  the  cosine  approximation  becomes  poor  about  30 
cm  away  from  the  mid-plane.  More  sources  are  needed  for  a 
smoother  axial  flux  profile  on  the  first  wall  of  a 
cylindrical  facility. 

A  similar  search  was  performed  for  five  point  sources  to 
determine  if  this  would  be  enough  for  a  smooth  flux  profile. 
Positions  and  Zg  were  added,  the  search  was  improved  as 
better  fits  were  obtained,  and  constraints  were  successively 
loosened.  The  constraints  were 

1 )  AZ  j »AZ 2 -AZ 3  *AZ 4  ,  S | — S 2  — S  3  BS 4  — S  5 , 

2)  AZ1*AZ4XAZ2=-‘AZ3  ,  S1=-Sa*S3=S4-S5  , 

3 )  AZ i =AZ 4  ^  AZ 2  * AZ  3 ,  S  j =S 5 2  *S  3  — S  4 , 

4)  AZ j -AZ 4 * AZ 2 “AZ 3 ,  S l =S5*S3*S3-S4 ,  and  finally 

5)  interactive,  making  slight  adjustments  in  any  source 
strength  or  position. 

The  result  of  this  search  is  shown  in  Figure  IV. 4;  the  flux 
profile  from  five  sources  is  obviously  a  very  close 
approximation  to  a  chopped-cosine  distributed  source. 
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The  RMS  error  for  this  search  was  less  than  0.1*  for  the 
central  80*  of  the  wall.  Except  for  the  energy  variation 
axially,  this  set-up  could  produce  outstanding  experimental 
results . 

Neutrons  from  an  axially  aligned  d— beam  would  produce 
flux  profiles  which  would  still  be  unsymmetric  in  energy, 
with  the  average  value  in  the  lower  end  of  the  blanket  about 
1.5  MeV  (~10*)  greater  than  in  the  upper  end.  The  effect  of 
this  variation  was  not  investigated,  but  it  would  tend  to 
cause  very  small  energy-related  differences  in  measurement 
results.  Another  detrimental  effect  of  using  few  sources  is 
ray-effects  in  neutron  transport  calculations. 34  The 
unattenuated  neutrons  striking  any  portion  of  the  blanket 
would  have  three  discrete  incoming  angles,  which  has  been 
shown  to  produce  erroneous  results  in  most  transport  codes. 
These  ray  effects  can  be  easily  treated  by  the  use  of  codes 
which  calculate  a  first— scatter  and  unscattered  neutron 
source  throughout  the  blanket.  Nevertheless,  the 
combination  of  ray  effects,  large  energy  transition  from  top 
to  bottom,  difficulty  of  alignment  of  three  or  more  targets 
with  a  charged  particle  beam,  unpredictable  transport  of  a 
charged  particle  beam  past  multiple  conducting  targets,  and 
necessary  bulky  cooling  structure  make  the  prospect  of  any 
"point  source"  combination  unattractive. 

Thus,  there  is  a  need  to  find  a  better  alternative.  The 
most  desireable  alternative  is  to  create  an  appropriate 


source  system  which  will  produce  a  continuous  flux 
distribution  on  the  first  wall.  The  solution  to  this  is 
addressed  in  the  following. 

Since  the  distance  from  the  axis  of  the  blanket  to  the 
first  wall  is  short,  and  due  to  the  long  mean— free— paths  of 
the  source  neutrons,  smoothing  of  the  neutron  flux  on  the 
first  wall  by  scattering  would  be  difficult  (but  not 
impossible)  in  the  proposed  FRBF.  Also,  too  much  scattering 
would  alter  the  flux  spectrum  so  it  would  not  represent  that 
incident  on  the  first  wall  of  a  real  fusion  reactor  blanket. 
The  solution  to  this  dilemma  is  to  create  an  actual  line 
source,  which  has  never  been  done  for  a  fusion  reactor 
blanket  experiment.  A  14  MeV  line  source  could  be  created 
by  using  an  axial  strip  of  TiT  and  a  modulated  deuteron 
beam.  The  d-beam  would  then  have  to  come  in  the  side  of  the 
cylindrical  blanket,  requiring  a  gap  in  the  structure  (as  in 
Ref.  24).  Studies  conducted  with  the  FBBF  have  shown  that 
the  neutron  flux  is  unaffected  by  a  discontinuity 
approximately  30  degrees  in  azimuth  away  from  measurement 
location.  An  experiment  to  validate  this  prediction  for  a 
vacuum  discontinuity  in  the  FBBF  will  be  performed  during 
the  unloading  of  the  current  blanket.  If  effects  away  from 
a  discontinuity  are  also  absent  from  the  FRBF  with  a  "slice" 
removed  from  the  side,  even  at  an  acceptably  greater 
azimuth,  a  greatly  improved  concept  for  fusion  reactor 
blanket  neutronics  evaluations  will  have  been  developed.  To 
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produce  this  line  source,  the  d-beam  will  be  directed  along 
the  target  by  a  modulating  electric  or  magnetic  field,  such 
that  it  produces  a  near  1 y  exact  chopped-cosine  flux 
distribution  in  the  center,  at  the  first  wall,  and  in  the 
blanket.  Figure  IV. 5  is  a  sketch  of  this  configuration, 
including  the  source  region,  blanket  components,  and  d-beam 
generator  system.  The  advantages  of  this  configuration  will 
be  many:  no  axial  or  radial  ray  effects,  only  very  gradual 
variation  of  energy  axially,  much  longer  lifetime  of  the 
target,  no  problem  with  axial  alignment  of  multiple  targets 
in  the  path  of  a  single  beam,  and  most  important,  a  source 
profile  which  matches  the  axial  flux  shape  in  the  blanket. 
However,  there  will  still  be  an  energy  and  yield  variation 
azimuthal ly;  their  effects  on  transport  predictions  require 
analysis.  Although  near  1 y-constant-density ,  long,  thin 
tritium  targets  can  be  produced,35  line  targets  could  be 
constructed  of  short  sections  of  TiT-plated  copper.  The 
tritium  density  of  these  targets  can  easily  be  measured  by 
detection  of  Bremsstrah 1 ung ,  which  will  ensure  a  consistent 
source  strength. 

The  neutron  source  scheme  described  above  requires 
evaluation  of  neutron  transport.  It  also  requires 
engineering  design  to  determine  d-beam  production  and 
control  methods,  source  strength  determination  methods  (a 
bank  of  recoil-alpha  detectors),  vacuum  system  requirements, 
and  beam  sweeping  and  control  requirements;  that  evalution 
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ketch  of  the  Fusion  Reactor  Blanket  Facility. 
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is  left  to  others.  The  presentation  of  the  requirements  for 
analysis  of  neutron  transport  in  this  proposed  fusion 
reactor  blanket  simulation  facility  follows. 


$ 


IV. 7.  Neutronics  Analyses 

The  blanket  design  and  source  arrangement  previously 
described  will  affect  the  transport  of  neutrons  into, 
through,  and  out  of  the  FRBF.  Conversely,  the  transport  of 
neutrons  in  the  FRBF  is  the  major  factor  defining  the 
selection  of  the  blanket  design  and  the  source  arrangement. 
Therefore,  the  critical  portion  of  the  research  reported  in 
this  paper  involves  neutronics  analysis  of  the  blanket/ 
source  configuration.  The  following  should  be  considered: 

•  transport  in  the  target  support  structure, 

•  effect  of  removing  part  of  the  azimuth, 

•  effect  of  azimuthal  energy  and  yield  variation. 

Transport  through  the  structure  surrounding  the  actual 
target  design  (backing,  coolant  channels,  support,  etc.) 
should  be  evaluated  to  determine  its  effect  on  the  neutron 
flux  and  spectrum  throughout  the  system.  Any  adverse  effect 
must  be  minimized  to  reduce  perturbations  which  cannot  be 
treated  exactly  in  2-dimensional  transport  codes.  However, 
the  construction  of  the  source  in  a  long  target  will  reduce 
the  thermal  flux  (per  unit  area)  by  about  two  orders  of 
magnitude,  thus  reducing  the  requirement  for  extensive 


cooling  structure.  The  almost  symmetric  vacuum  tube  located 
a  distance  from  the  source  also  eliminates  some  of  the 
asymmetry  of  previous  sources  caused  by  support  and  cooling 
structure.  Additionally,  the  target  can  be  constructed  so 
that  it  is  convex,  thus  offering  no  long,  straight  path  for 
neutron  attenuation  in  the  radial  direction.  Then, 
absorption  and  scattering  of  the  source  neutrons  will  be 
almost  negligible,  and  will  not  contribute  significantly  to 
asymmetry  of  flux  distributions  in  the  blanket.  Therefore, 
the  source  region  can  be  treated  as  a  symmetric  target,  and 
extensive  neutronics  analyses  of  the  target  system  are  not 
required  for  this  research,  although  these  analyses  must  be 
done  for  the  evaluation  of  experiments  in  a  final  system. 

The  effect  of  removing  part  of  the  blanket  for  the  d— beam 
port  must  be  analyzed  to  determine  if  it  severely  affects 
the  flux  or  spectrum  many  mean— f ree-paths  from  the 
disconcinuity .  An  azimuthal  segment  without  large  flux 
deviations  spanning  only  45°  to  135°  from  the  asymmetry  will 
allow  an  experimental  arc  of  about  40  to  100  cm  in  the 
breeding  zone  and  150  cm  at  the  outer  edge  of  the  shield. 

To  reduce  the  costs  of  multiple  evaluations,  the  effect  of 
the  asymmetry  also  must  be  minimized  to  maintain 
approximately  1-dimensional  geometry.  Finally,  the  effect 
of  the  azimuthal  asymmetry  of  flux  and  spectrum,  due  to  the 
anisotropic  neutron  source,  must  be  determined. 
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The  previously  mentioned  effects  were  determined  by 
comparing  1-D  transport  predictions  (symmetric  blanket  with 


isotropic  source)  to  2-D  predictions  (asymmetric  blanket 
with  anisotropic  source).  The  method  used  to  make  those 
predictions  and  comparisons  is  presented  in  the  following 
chapter . 


V.  NEUTRON ICS  ANALYSES  OF  TWO  FRBF  MODELS 


The  preceding  chapters  explained  the  requirements  for 
more  and  better  fusion  blanket  experiments,  and  described  a 
proposed  new  facility,  the  FRBF.  for  conducting  those 
experiments.  To  insure  a  new  experiment  or  facility  will 
yield  improved  data,  it  should  be  analyzed  to  help  eliminate 
design  deficiencies.  This  chapter  describes  the  analysis  of 
neutron  transport  through  the  two  models  which  were  selected 
for  analyzing  the  FRBF  design.  To  perform  such  an  analysis, 
a  neutron  transport  code  is  required.  The  code  requires  a 
correct  source  description  for  input,  an  appropriate  set  of 
group  constants  to  describe  the  materials’  transport 
characteristics,  and  correct  modeling  of  the  system 
geometries.  The  neutron  transport  code  used  for  this 
analysis  is  discussed  first,  followed  by  a  description  of 
the  method  used  to  generate  the  isotropic  and  anisotropic 
source  files  for  input  to  the  code.  Then,  the  source  and 
production  of  the  group  constants  used  in  this  analysis  is 
discussed,  along  with  a  presentation  of  the  importance  of 
inelastic  scattering  reactions  in  fusion  materials. 

Finally,  the  use  of  the  different  models  for  analyzing 
transport  through  the  FRBF  is  described.  The  order  of 


investigation  of  the  transport  is  presented  with  a 
discussion  of  the  output  data,  which  was  used  to  produce  the 
information  presented  in  Chapter  VI,  Neutronics  Results. 

• 

V.l.  DORT ,  A  Transport  Code 

A  new  discrete  ordinates  theory  (D.O.T.)  code  was  used  to 
analyze  effects  on  neutron  transport  predictions  due  to  the 
proposed  geometry  and  source  configuration  of  the  FRBF.  The 
code,  DORT  1.0,  was  produced  by  the  Reactor  Shielding 
Information  Center  (RSIC)  at  ORNL.  This  code  is  an  upgraded 
version  of  DOT  5.1,  which  was  the  next  generation  of  the  DOT 
series.36  The  previous  version,  DOT  4.3,  is  in  wide  use  at 
many  laboratories  and  industries,  and  was  used  by  the  BCSS 
for  blanket  design  and  selection  studies  and  by  RSIC  for 
deep  penetration  (shielding)  problems.  Prior  to  this 
research,  no  version  of  DOT  was  available  at  Purdue 
University.  RSIC  gave  this  DORT  package  to  Purdue  for 
modification  and  validation;  it  was  converted  for  this 
research  to  run  on  the  Cyber  205  vector  supercomputer,  and 
it  was  thoroughly  tested  with  documented  sample  problems 
prior  to  use  for  FRBF  analyses.  The  sample  problems 
included  those  supplied  with  the  DOT  4.3  and  DORT  code 
packages,  and  two  problems  using  EPRI-Cell  cross  sections. 

In  most  cases,  the  Purdue  University  version  of  DORT 
produced  results  which  were  identical  (effective  k  and  flux 
values)  to  the  reference  (Cray)  cases  from  RSIC.  Using  the 
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EPRI-Cell  benchmark  cross  sections,  DORT  produced  almost 
identical  results  (differences  ~0.01SS)  compared  to  TWODANT87 
(D.O.T.)  and  2DB38  (diffusion  theory)  results  computed  at 
Purdue . 

The  DORT  code  is  a  multi-capability  program  offering  many 
alternatives  for  input,  output,  source  types,  mesh  spacing, 
quadrature  sets,  etc.  It  allows  variable  mesh  spacing  and 
quadrature,  and  accepts  a  variety  of  sources  (alone  or 
combined);  internal  and/or  external  boundary  sources, 
distributed  source,  first-collision  source,  etc.  An 
additional  feature  is  a  capability  for  removing  negative 
fluxes  which  result  from  truncation  of  a  polynomial 
expansion  of  the  scattering  terra  in  the  D.O.T.  treatment. 
This  function  is  especially  useful  where  there  are  large 
void  regions.  Another  extremely  useful  feature  is  the 
capability  to  save  flux,  flux  moment,  and  boundary  source 
files;  then  to  use  those  as  input  or  re-start  data  for 
subsequent  problems.  This  feature  was  used  for 
investigations  of  R-Z  and  R-0  geometries,  with  either  a 
distributed  or  a  boundary  source.  The  generation  of  these 
sources  is  explained  in  the  following  section. 


V.2.  The  Anisotropic  Neutron  Source 

Production  of  neutrons  by  the  T(d,n)a  reaction  was 
described  earlier  as  being  dependent  on  the  energy  of  the 
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deuteron  beam.  Slowing  down  of  the  beam  ions  in  the  target 
causes  a  distribution  of  neutron  energies  at  any  angle  of 
emission.  The  treatment  of  this  was  described  clearly  in 
Ref.  32,  which  includes  tables  of  the  energy-angle 
correlation.  The  Fortran  program  included  in  Ref.  32  was 
modified  for  this  research,  and  is  listed  in  Appendix  B.  It 
employs  the  slowing  down  and  energy— angle  correlation  to 
produce  a  table  of  neutron  emission  probabilities  as  a 
function  of  energy  group  and  angle  interval.  However,  the 
angular  treatment  was  improved  for  a  better  approximation  to 
the  integral  over  emission  angle  in  the  forward  and  backward 
directions  (near  0°  and  180°).  Table  7.1  is  one  result  of 
this  program,  with  the  distribution  computed  for  a  beam 
energy  of  150  KeV.  Note  that  the  isotropic  probability 
between  0  and  20°  would  be  0.03015,  while  the  anisotropic 
probability  is  0.03155  (~5%  difference).  This  program  also 
predicts  the  theoretical  emission  rate  of  neutrons  per 
milliamp  deuteron-beam  current;  Figure  V.l  is  an 
illustration  of  this  relationship.  These  values  are  in  good 
agreement  with  others,39  as  are  the  probabilities  versus 
angle  and  energy,  which  were  used  to  generate  sources  for 
the  neutronics  evaluations. 

The  results  in  Table  V.l  were  used  in  two  ways  for  input 
to  the  transport  calculations.  For  1-D  predictions,  the 
probabilities  for  each  energy  group  were  summed  over  all 
directions,  so  the  input  spectrum  (at  the  bottom  of  Table 
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♦Read  1 . 000e-2  as  1.000- 10_‘°,  values  are  probability  of  emission  of  a  d-T 
neutron  into  the  solid-angle  interval  and  energy  group.  Sums  of  the 
columns  gives  the  spectrum,  while  sum  of  a  row  gives  the  total  probability 
of  emission  of  a  neutron  into  a  solid-angle  interval. 
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Figure  V.l.  Yield  of  d-T  Neutrons  vs  d-beam  Current.  Based 
on  numerical  integration  of  the  slowing  down  equation  for 
deuterons  in  TiT2 .  assuming  the  TiT  layer  is  thick  enough  to 
stop  the  deuterons. 

V.l  was  an  isotropic  distribution  of  those  five  groups.  For 
the  2-D,  asymmetric  blanket  analyses,  a  second  Fortran 
program  was  written  which  correlated  the  direction  cosine  of 
the  face  of  a  mesh  space  (cos <P)  and  the  direction  cosines  of 
the  quadrature  set  (m  and  7?)  with  the  angle  of  emission  of 
the  neutron  from  the  target  (0).  The  equation  which  relates 
these  angles  and  directions  is 

cos  0  =  M'  cos^  -  77‘sin<£.  (5.1) 

The  direction  cosine  with  respect  to  the  radial  direction  is 
fi ,  and  the  direction  cosine  with  respect  to  the  azimuthal 
direction  is  77.  The  program  selects  from  the  table  a 
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obtained  for  the  transport  computations.  The  group  constant 
library  used  to  begin  the  transport  investigation  was 
produced  by  R.  Roussin  at  RSIC  from  the  VITAMIN-E  (ENDF/B-V) 
library.  He  used  a  series  of  programs  to  select  the  fifteen 
isotopes  and  to  collapse  the  174-neutron-group  data  base  to 
50  groups,  with  a  P3  scattering  expansion.  Upscatter  was 
not  included,  and  no  neutron  production  cross  sections  were 
included.  Tritium  production  cross  sections  were  also  not 
provided  with  this  data  set.  The  energy  bounds  and  lethargy 
widths  are  listed  in  Appendix  A.  The  collapsed  (but 
unmixed)  library  was  transmitted  on  a  tape  for  use  at  Purdue 
University.  A  second  library,  including  the  same  isotopes, 
but  also  including  tritium  breeding  constants,  was  generated 
by  I.  Dilber  from  the  MARS/VI TAM I N-E  package  recently 
activated  on  the  Purdue  University  computer  system.  The 
final  transport  calculations  were  completed  using  this 
forty— one— group  library,  after  checking  for  consistency  with 
the  RSIC— supp 1 ied  constants.  However,  the  constants  from 
either  source  required  mixing  and  manipulation  prior  to 
input  to  DORT. 

The  cross  sections  were  prepared  for  input  to  the  DORT 
code  by  processing  them  through  GIP,  another  program  in  the 
DORT  package  of  codes.  GIP  accepts  isotope-ordered  cross 
sections,  combines  them  in  accordance  with  an  input  mixing 
table,  and  its  output  is  a  group-ordered  library  of  mixed 
constants  suitable  for  use  by  DORT,  ANISN,  DOT  3,  or  DOT  4. 


56 


GIP  was  used  to  prepare  the  material  constants  by  using  the 
number  densities  listed  in  Appendix  A.  Additionally,  pure 
8Li  and  7Li  constants  were  supplied  for  Blanket  HeC  to 
compute  reaction  rates  for  tritium  breeding  comparisons.  In 
Blanket  LLC,  because  T7  contributes  only  0.2%  to  the 
breeding,1  tritium  breeding  was  not  computed  for  the 
individual  isotopes,  only  for  the  total  material.  The 
macroscopic  constants  which  were  produced  by  GIP  were  then 
used  for  the  transport  predictions  with  the  models  described 
in  the  following  section. 

This  is  an  appropriate  place  to  emphasize  the  importance 
of  the  inelastic  reaction  to  neutronics  in  fusion  reactor 
materials.  Table  V.2  lists  elastic  and  inelastic  scattering 
cross  sections  for  a  few  important  isotopes  in  fusion 
reactor  blankets,  shields,  first  walls,  and  magnets.  The 
two  types  of  reactions  occur  on  the  same  order  of  magnitude. 
However,  the  energy  loss  in  the  two  reactions  is  very 
different.  Figure  V.2  illustrates  the  energy  loss  due  to 
inelastic  reactions  for  lead  and  copper.  These  cross 
sections  were  taken  from  LIB-IV.41  On  the  average,  in  a 
single  collision,  a  14  MeV  neutron  loses  12.1  MeV  in  Pb 
(86%);  in  contrast,  only  0.2%  is  lost  in  an  elastic 
collision.  The  energy  which  is  lost  in  these  inelastic 
collisions  is  deposited  locally  due  to  short-range  gamma- 
rays  emitted  in  the  decay  of  excited  nuclei.  Thus,  a  large 
part  of  the  power  density  in  fusion  reactor  materials  is  due 
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Table  V.2.  Scattering  Cross  Sections  for  Fusion  Reactor 
Materials 


Isotope 

Energy 

(MeV) 

Elastic 

Scatter 

Inelastic 

Scatter 

Energy 

Left 

•Li 

14 

0.97* 

0.372 

0.936 

10 

1.25 

0.589 

0.888 

1 

1.05 

0.000 

0.769 

7Li 

14 

1.01 

0.370 

0.929 

10 

1.41 

0.538 

0.885 

1 

1.38 

0.  179 

0.763 

Fe 

14 

1.23 

0.782 

0.992 

10 

1.98 

1.47 

0.993 

1 

2.  13 

0.325 

0.997 

Cu 

14 

1.39 

0.894 

0.995 

10 

2.  10 

1.56 

0.994 

1 

2.92 

0.269 

0.974 

Pb 

14 

2.91 

0.359 

0.998 

10 

3.36 

2.32 

0.997 

1 

4.82 

0.223 

0.992 

♦Cross  section  in 

barns.  Energy  left  is  the  average 

of  the  neutron's 

original  energy  which  remains  after 

i  one  elastic  collision.  This  can  be  contrasted  to  I 

the  much 

lesser  fraction  remaining  (0.05  to 

0.20) 

after  one  inelastic  collision. 
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lastic  Scattering  Cross  Sections  of  Copper 
MeV  Neutrons.41  This  figure  demonstrates 
losses  due  to  inelastic  scattering  of  14 
lead  and  copper,  and  of  6—10  MeV  neutrons  j 
rage  energy  of  the  neutron  following  the 
ated.  Thus,  the  average  14MeV  neutron 
ticaily  in  lead  deposits  12.1  MeV  (86%)  of 
e  lead  nucleus.  This  energy  is  dissipated 
products . 
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to  the  inelastic  scattering  process.  Another  demonstration 
of  the  importance  of  this  reaction  is  shown  in  Figure  V.3, 
which  compares  the  energy  loss  due  to  elastic  scattering 
with  the  loss  due  to  inelastic  scattering  in  one  FRBF  model. 
This  figure  represents  the  product  of  the  isotopic  reaction 
rate  with  the  average  energy  loss  rate  for  the  inelastic  or 
elastic  reaction  (based  on  results  from  Chapter  VI).  It 
demonstrates  that  the  poswer  density  due  to  inelastic 
scattering  will  be  a  factor  of  ten  greater  than  that  due  to 
elastic  scattering  in  a  fusion  reactor. 

The  important  conclusion  to  be  drawn  from  this  is  that 
knowledge  of  inelastic  cross  sections  is  critical  for 
predicting  the  location  and  density  of  energy  deposition  in 
a  fusion  reactor  system.  However,  accurate  determination  of 
these  cross  sections  is  very  difficult  due  to  the  high- 
energy  sources  required,  the  multitude  of  competing 
reactions,  and  the  spread  of  energy  and  direction  of 
neutrons  emitted  in  the  various  reactions.  Therefore,  these 
complicated  scattering  kernals  or  two-dimensional  matrices 
of  inelastic  reactions  are  not  accurately  determined. 

Neither  can  the  power  distribution  in  a  fusion  reactor  be 
accurately  predicted,  based  on  available  cross  section 
information  and  current  computational  methods.  As  this 
knowledge  is  critical  for  accurate  prediction  of  neutronics 
and  heat  production  for  fusion  reactor  design,  experimental 
studies  are  required  to  increase  confidence  in  the  data  and 
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capability  to  accurately  compute  this  information. 


V.4.  Computational  FRBF  Models 

Four  different  models  were  required  for  the  evaluations 
for  this  research.  First,  an  R-Z,  2-D  calculation  was  done 
using  a  P3  expansion  of  the  scattering  cross  sections  and  an 
S16  quadrature  set.  The  results  of  this  calculation  were 
used  to  estimate  DB2  leakage  for  input  to  the  R-0 
computations,  and  for  demonstrating  the  benefits  of  proper 
source  selection  (shape).  The  blanket  was  divided  on  the 
plane  and  axis  of  symmetry;  it  was  treated  only  from  the 
mid-plane  up  and  the  axis  out,  with  axial  mesh  spacing  every 
5  cm  (about  1/2  mean-f ree-path  for  fusion-energy  neutrons). 
The  bottom  and  left  (center)  boundary  conditions  were 
reflective,  while  the  top  and  right  (outer)  boundary 
conditions  were  void. 

The  radial  geometries  of  the  R-Z  and  R— 8  computations 
were  identical;  the  spacing  was  given  in  Figures  IV. 2. A  and 
B.  The  breeder  zones  were  divided  into  5  cm  spacing  (again 
''■'1/2  MFP) ,  while  the  reflector  or  plenum  and  the  shield  were 
of  slightly  coarser  mesh.  The  azimuthal  intervals  were  10°, 
with  only  one-half  azimuth  of  the  blanket  modeled  (18 
intervals).  The  boundary  conditions  for  the  azimuthal 
surfaces  were  reflective,  the  left  (center)  boundary 
condition  was  cylindrical,  and  the  right  (outer)  boundary 


condition  was  void. 
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The  DB2  leakage  from  the  R-Z  computation  was  then  used  in 
a  1-D  calculation  to  determine  the  flux  and  tritium 
production  for  a  symmetric  system  with  an  isotropic 
distributed  source.  These  results  could  then  be  compared  to 
the  results  from  a  computation  with  an  asymmetric  model. 

The  model  for  this  2-D  computation  included  a  10°  segment 
which  contained  a  vacuum  (air  at  10"3  atm).  Also  included 
was  a  SS  wall  in  the  next  10°  segment  between  the  vacuum 
tube  and  the  first  wall.  The  vacuum  region  simulated  the 
section  removed  from  the  blanket  for  the  incoming  d-beam, 
and  the  SS  wall  represented  the  containment  necessary  to 
maintain  a  vacuum  for  the  neutron  generating  system.  The 
radially  outermost  mesh  space  of  the  first  azimuthal 
interval  contained  shield  material,  rather  than  vacuum,  to 
represent  the  closure  of  the  vacuum  system  and  shielding  and 
instrumentation  which  will  be  required  for  the  neutron 
generating  system.  Figures  V.4.A  and  B  illustrate  the 
radial  zone  boundaries  and  the  asymmetry  for  the  blanket 
models.  The  10°  intervals  represent  a  SS  wall  3.8  cm  thick 
at  the  first  wall,  a  7.7  cm  vacuum  gap  (2  times  the  half¬ 
geometry)  at  the  first  wall,  and  a  35  cm  section  at  the 
outer  radius  of  the  shield.  These  large  dimensions  can  be 
considered  upper  limits  on  the  expected  assymetries;  they 
should  produce  maximum  deviations  in  the  results. 
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Two  more  models  were  used  to  estimate  the  differences 
between  experimentally  validated  blanket  physics  and  actual 
fusion  reactor  physics.  The  first  was  similar  to  the  1-D 
model  discussed  previously;  however,  the  axial  leakage  was 
zero  (as  in  an  infinite  system  or  the  center  of  a  long 
system).  The  second  model  was  also  1-D  with  zero  leakage, 
but  instead  of  a  source  distributed  in  a  thin  target  in  the 
center,  a  plasma  source  was  distributed  in  a  16-cm-radius 
vacuum  region.  These  two  models  were  used  fcr  both  blanket 
compositions  to  determine  flux  profiles  and  tritium 
production  rates.  These  four  models  were  used  to  predict 
neutron  transport  through  both  Blankets  HeC  and  LLC ;  the 
results  are  displayed  in  the  following  chapter. 


VI.  NEUTRONICS  RESULTS 


The  preceding  chapters  discussed  the  requirement  for  a 
new  facility  for  testing  neutron  transport  in  fusion  reactor 
blankets,  a  conceptual  design  for  a  facility  which  will  meet 
that  need,  and  the  theory  behind  the  design  of  that 
facility — the  FRBF .  Chapter  V  described  a  neutronics  code 
to  analyze  that  proposed  design,  the  models  and  sources 
necessary  for  the  analysis,  and  the  group  constants  which 
could  be  used  to  predict  the  transport  of  neutrons.  This 
chapter  presents  and  compares  the  results  of  that  neutronics 
evaluation.  The  results  are  presented  in  graphical  and 
tabular  forms  as  flux  calculations  and  tritium  breeding 
results  from  the  various  computations.  This  chapter 
discusses  the  neutronics  results  in  the  same  order  as  the 
calculations  were  discussed  in  section  V.5:  R-Z,  1-D  R-0 
versus  2-D  R -6  with  asymmetry,  and  1-D  R-0  central  source 
versus  1-D  R-0  plasma  source.  The  results  of  each  model  are 
presented  as  the  appropriate  flux  comparison,  tritium 
breeding  comparison,  etc.  All  results  have  been  normalized 
to  a  per-source-neutron  basis  for  each  model.  Each  result 
for  Blanket  HeC  is  followed  immediately  by  the  result  for 
Blanket  LLC  in  each  section.  An  attempt  is  made  to  present 


the  results  in  the  most  informative  manner  for  the  aspect  to 
be  demonstrated,  as  group-flux  profiles,  position-dependent 
spectra,  flux  ratios,  or  tritium  production  profiles.  Since 
the  R-Z  analysis  was  required  for  buckling  corrections  for 
input  to  the  R-6  predictions,  and  to  demonstrate  the 
adequacy  of  the  cosine-shaped  axial  source,  it  is  presented 
first . 

VI. 1.  Results  of  2-D,  R-Z  Analyses 

The  R— Z  computations  were  completed  to  satisfy  three 
separate  needs: 

•  buckling  corrections  were  needed  for  the  R-9 
geometries,  to  obtain  more  realistic  transport  results; 

•  axial  source  and  group-flux  profiles  were  needed  to 
demonstrate  the  suitability  of  the  cosine-shaped  line 
source;  and 

•  leakage  fluxes  were  needed  to  calculate  dose  rates 
in  the  FRBF  and  surrounding  areas. 

The  buckling  coefficients  which  were  obtained  for  Blankets 
HeC  and  LLC  were  used  as  input  to  the  R-0  predictions.  No 
presentation  of  those  values  is  given  here,  as  they  would 
demonstrate  nothing  significant.  However,  the  axial  flux 
profiles,  which  develop  due  to  a  chopped— cosine— shaped  line 
source,  are  extremely  important,  as  this  is  the  basis  of  one 
of  the  major  advantages  of  the  FRBF.  Normalized  axial  flux 
profiles,  which  should  approximately  match  the  source 
profile,  are  shown  in  Figures  VI. 1. A  and  B.  These  figures 


RELATIVE  FLUX 


DISTANCE  FROM  MID-PLANE  (CM) 


Figure  VI. 1. A.  Axial  Group-flux  Profiles  in  Blanket  HeC. 
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also  include  the  source  profile,  which  is  a  chopped-cosine 
function  in  the  axial  dimension  (same  for  all  source 
groups).  The  approximate  equality  of  the  flux  profiles  with 
the  source  profile  demonstrates  the  suitability  of  the 
cosine-shaped  line  source  to  this  geometry  (Figure  IV. 3 
showed  how  a  point  source,  which  is  a  delta  function  at  the 
mid-plane,  transitions  to  a  gaussian  at  the  first  wall;  the 
flux  profile  does  not  match  the  source  profile).  The 
variation  of  flux  profile  from  the  cosine  distribution  at 
the  top  (and  bottom)  of  blanket  HeC,  near  the  shield,  is 
from  isotropic,  elastic  scattering  of  lowei — energy  neutrons 
in  the  top  shield.  This  strongly  supports  the  previously 
stated  belief  that  spectrum  transition  will  be  a  large 
factor  in  a  fusion  blanket.  Figures  VI. 2. A  and  B  show  the 
spectra  in  various  regions  of  the  blankets,  and  Figures 
VI. 3. A  and  B  present  the  radial  flux  profiles  for  groups  3 
(~14  MeV),  7  (~10  MeV),  11  (~1  MeV) ,  and  14  (~300  KeV) . 

These  two  figures  also  help  demonstrate  the  flux  transition 
within  the  blanket.  The  neutron  spectrum  gets  'softer'  as 
the  flux  travels  through  the  blanket  (due  to 
downscattering).  However,  Figures  VI. 2. A  and  B  demonstrate 
that  the  spectral  shapes  are  almost  identical  axially;  there 
is  no  axial  energy  dependence  of  the  flux.  Finally,  the 
axial  tritium  production  predicted  for  the  blankets  is 
illustrated  in  Figure  VI. 4. 
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Figure  VI. 2. A.  Neutron  Spectra  in  Blanket  HeC ,  R-Z. 
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Figure  VI. 2. B.  Neutron  Spectra  in  Blanket  LLC ,  R-Z 
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Figure  VI. 3. A.  Radial  Group-flux  Profiles  in  Blanket  HeC 
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VI. 2.  Comparison  of  1-D  and  2-D  R-8  Results 

To  determine  the  effects  of  the  asymmetry  of  this  design 
of  the  FRBF ,  a  2-D  transport  computation  was  performed. 

This  section  presents  the  results  of  the  comparison  of  the 
2-D  asymmetric  models  to  the  1-D  symmetric  models,  both  with 
isotropic  distributed  sources.  The  first  figures,  Figures 
VI. 5. A  and  B,  present  the  flux  ratios  for  various  groups  at 
mid-breeder  as  they  vary  with  azimuth.  The  flux  ratio  is 
defined  as  the  value  of  flux  obtained  from  the  1-D  model 
divided  by  the  value  from  the  2-D  model.  These  figures 
demonstrate  that  in  an  arc  of  about  90°  to  180°  away  from 
the  asymmetry,  the  high-energy  2-D  flux  is  within  ~Z%  of  the 
1-D  prediction.  The  lowei — energy  ratios  deviate  more  from  1 
in  the  Blanket  LLC ,  which  is  an  unexpected  result.  This  may 
be  due  to  a  lack  of  capture  of  neutrons  in  the  void  region, 
to  scattering  of  excess  high-energy  flux  in  the  outer 
boundary  of  the  void,  to  scattering  in  the  vacuum  wall,  or 
some  other  effect.  It  also  may  be  due  to  truncating  the 
expansion  of  the  scattering  matrix  at  P3 ,  or  to  the  negative 
source  removal  in  DORT.  This  deviation  will  require  further 
study  to  determine  its  cause  and  change  the  design  to 
ameliorate  it  (see  section  VIII. 2,  Recommendations).  Radial 
flux-ratio  profiles  are  then  shown  in  Figures  VI. 6. A  and  B 
for  various  groups  at  0“1O5°,  where  the  2-D  results  are 
close  to  the  1-D  results.  These  figures  demonstrate  that 
the  relationship  between  the  two  predictions  becomes  closer 
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to  1  as  the  neutron  flux  travels  deeper  into  the  blanket. 
The  presence  of  the  void  region  thus  has  a  strong  effect  on 
the  neutron  source  in  the  inner  portion  of  the  blanket. 
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Figures  VI. 7. A  and  B  present  the  tritium  production  rates 
for  various  azimuthal  segments  of  the  blankets.  Since 
tritium  from  *Li(n,T)a  reactions  occurs  at  lower  energy,  the 
effect  of  the  deviation  of  the  lower — energy  neutron  flux  due 
to  the  asymmetry  is  clear  in  Figure  VI. 7. B.  Figure  VI. 7. A 
shows  that  T7  is  almost  unaffected  by  the  asymmetry,  because 
the  high-energy  flux  which  produces  T  in  the  7 Li (n , n ’T)a 
reaction  is  not  affected  by  the  asymmetry.  This  is  because 
the  highest-energy  neutrons  travel  radially  outward  until 
they  are  elastically  scattered,  ine 1 ast ica 1 1 y  scattered,  or 
absorbed.  Only  the  elastically  scattered  neutrons  stay  in 
the  high  energy  groups,  and  the  average  angle  of  deflection 
in  high  energy  is  near  zero.  The  average  cosine  of  the 
scattering  angle  for  iron,  for  instance,  is  0.83  at  10  MeV. 
Thus,  high-energy  neutrons  would  have  a  very  small  chance  of 
scattering  90°  around  the  blanket,  from  the  void  region. 

The  previous  comparisons  do  not  include  the  effects  of  the 
anisotropic  source,  which  does  have  a  slight  effect  on  the 
flux  distribution  in  the  blanket. 

2-D  calculations  with  anisotropic  sources  demonstrated 
that  the  flux  in  any  region  generally  varied  as  the  source 
in  groups  2  and  3  varied.  However,  these  computations  were 
not  completed  in  a  consistent  manner,  due  to  an 
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incompatibility  between  the  DORT  treatment  of  the  quadrature 
and  boundary  flux  input  file  in  DORT.  They  did  show  that  if 
the  flux  in  group  3  was  greater  at  an  azimuth,  the  flux  in 
the  lower  groups  would  also  be  greater  (but  not  as  strongly) 
at  that  azimuth.  Thus,  since  the  source  strength  varies  by 
about  ± 2 X  around  90°  from  the  d-beam,  the  variation  in  flux 
due  to  anisotropy  will  be  less  than  ±2%  in  this  region  (60°— 
120°).  Because  this  effect  can  be  accounted  for  and  is  not 
affected  by  the  assymetry,  the  2-D  comparisons  just 
presented  are  adequate  to  demonstrate  the  range  of 
deviations  due  to  the  design  of  the  FRBF . 

VI. 3.  Central  Source  vs  Plasma  Source  Results 

Neutronics  evaluations  of  the  FRBF  conceptual  design  with 
a  plasma  source  were  conducted  to  determine  the  differences 
between  predictions  based  on  FRBF  measurements  and  actual 
conditions  in  a  fusion  reactor  (plasma)  test.  Figures 
VI. 8. A  and  B  show  the  neutron  spectra  in  the  first  wall  and 
mid— breeder  for  the  central  (FRBF)  source.  Differences 
between  these  predictions  and  predictions  with  a  plasma 
source  are  shown  in  Figures  VI. 9. A  and  B.  The  large 
deviation  at  about  14  MeV  is  due  to  the  different  source 
spectra  in  the  two  models.  The  deviation  in  the  lower 
energies  is  not  due  to  this  spectral  difference  though, 
since  the  cross  sections  for  the  higher  energy  groups  are 
nearly  constant.  Also,  the  previous  discussion  demonstrated 


LETHARGY  FLUX  ( N/CM2-S0URCE  NEUTRON) 


LETHRR6Y  FLUX  ( N/CM2-SQURCE  NEUTRON) 


PERCENT  DIFFERENCE 


89 


that  the  inelastieal ly  scattered  neutrons  would  have 
approximately  the  same  energy,  independent  of  the  higher 
energy  group.  The  deviation  in  neutron  flux  in  the  blanket, 
then,  can  only  be  caused  by  a  difference  in  angular 
distribution  of  source  neutrons  with  respect  to  the  radial 
direction.  The  angular  flux  at  the  first  wall  due  to  the 
central  source  is  very  different  from  the  fusion  source, 
because  neutrons  from  the  central  source  must  travel 
radially  (perpendicular  to  the  azimuthal  direction).  In 
contrast,  neutrons  leaving  the  plasma  from  the  outer  edge 
can  exit  tangentially  bo  the  plasma,  thus  striking  the  first 
wall  at  an  angle  far  from  perpendicular. 

Figure  VI. 10. A  and  B  next  show  the  predicted  radial 
tritium  production  profiles:  again,  the  difference  is  due  to 
the  longer  path  length  of  some  of  the  high-energy  neutrons 
in  the  plasma  model.  Loss  of  neutrons  by  exponential 
attenuation  due  to  inelastic  scattering  is  significantly 
different  in  the  two  models,  as  is  the  production  of  lowei — 
energy  neutrons  from  this  reaction.  These  comparisons 
clearly  demonstrate  that  the  neutronics  conditions  in  a  test 
system,  even  with  an  axially-distributed  central  source, 
cannot  reproduce  the  conditions  in  a  reactor  (of  any  size) 
which  contains  a  plasma.  Thus,  for  any  confidence  in 
predicting  tritium  breeding,  energy  deposition,  power 
density,  and  material  damage  in  a  fusion  reactor  design 
prior  to  actual  testing;  the  neutronics  methods  must  be 
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Figure  VI. 10. A.  Central  vs  Plasma  Source:  Radial  Tritium 
Production  in  Blanket  HeC. 


TRITIUM  PRODUCTION  ( T/SOURCE  NEUTR0N/CM2 


validated  and  the  accuracy  of  data  libraries  must  be 


confirmed. 


The  FRBF  could  be  used  for  this  validation  of  the 
neutronics  methods  and  data  libraries.  The  next  chapter 
discusses  the  advantages  of  constructing  an  FRBF,  including 
a  discussion  of  constructing  the  facility  at  Purdue 
University.  It  also  compares  in  several  ways  the  FRBF  with 
past  and  planned  fusion  reactor  neutronics  experiments. 


VII.  ADVANTAGES  OF  THE  FRBF 


Results  of  the  investigations  of  neutronics  in  the 
proposed  FRBF  design  demonstrate  that  the  facility  can  be 
used  to  perform  accurate  neutronics  experiments,  with  only 
slight  deviations  due  to  asymmetry.  If  these  deviations  are 
predicted  and  considered  when  examining  experimental 
results,  they  should  not  contribute  significantly  to 
deviations  from  neutronics  computations.  Thus,  the 
construction  of  this  facility  should  allow  accurate 
representation  in  2-D  analyses  for  neutronics  predictions  of 
power  distribution  and  tritium  breeding.  The  major,  very 
important,  advantage  is  that  3-D  analyses  would  not  be 
required,  because  the  design  of  the  facility  allows 
separation  of  the  axial  flux  dependence  in  space  and  energy. 
This  chapter  describes  some  of  the  other  advantages  which 
the  FRBF  would  have,  compared  to  past  and  planned 
experiments.  It  also  addresses  additional  advantages  of 
constructing  the  facility  at  Purdue,  and  it  describes  other 
issues  which  should  be  considered.  One  of  the  major 
advantages  of  the  FRBF  would  be  flexibility  for  multiple 
experiments,  and  the  low  cost  which  this  flexibility  makes 


One  of  the  problems  with  past  FRB  experiments  has  been 
that  the  blanket  assemblies  were  one-of-a-kind,  or  the 
blankets  which  were  designed  for  flexibility  had  geometries 
which  prevented  accurate  modeling  (refer  to  Chapter  II). 

The  FRBF  could  eliminate  these  problems  because  it  was 
conceived  with  flexibility  as  a  design  goal,  as  was  the 
FBBF.  By  being  constructed  of  concentric  walls  and  an 
azimuthal ly  and  radially  segmented  shield,  it  offers  the 
opportunity  for  varying  configurations,  materials, 
densities,  etc.  axially,  radially,  and/or  azimuthal ly.  The 
previously  discussed  results,  for  blankets  having  different 
compositions  of  similar  isotopes,  demonstrated  that  spectra 
and  flux  values  were  similar  for  the  two  models  studied. 

This  makes  possible  a  further  flexibility  for  experiments: 
it  appears  the  two  ha  1 f s  of  the  FRBF  could  be  constructed  of 
slightly  different  materials  or  configurations,  and  neutron 
transport,  power  distributions,  and  tritium  breeding  could 
be  studied  in  each  at  about  60-120  degrees  from  the  d-beam 
direction  (the  strong  azimuthal  variation  seen  in  the  flux 
near  the  void  and  SS  wall  would  be  much  less  pronounced  near 
the  interface  between  two  similar,  but  different,  blankets). 
Additionally,  a  segmented  shield  would  permit  removal  of  a 
section  for  installation  and  testing  of  a  real  fusion 
blanket  module,  such  as  the  Lithium  Blanket  Module.30  The 
composition  of  the  rest  of  the  FRBF  blanket  would  reduce 
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flux  variations  at  the  boundaries  of  the  test  module.  This 
designed-in  flexibility  is  one  of  the  reasons  the  FRBF  could 
be  used  to  conduct  multiple  tests  at  reduced  cost. 

The  financial  advantage  of  this  flexible  design  is  easily 
demonstrated.  The  original  cost  of  the  FBBF  project, 
including  design,  fabrication  and  construction,  and 
conducting  the  first  experiment,  was  about  1/2  million 
dollars.42  In  contrast,  the  cost  of  the  current  project  to 
re-design  and  replace  the  breeding  blanket  is  about  25 
thousand  dollars;  a  factor  of  20  less.43  The  projected  cost 
for  the  most  complicated  blanket  design  is  not  expected  to 
exceed. 50  thousand  dollars,  or  1/10  the  cost  of  the  original 
facility.  These  costs  do  not  include  the  cost  of  fuel  for 
the  facility,  which  is  in  the  form  of  U02  fuel  rods  obtained 
on  a  loan  basis.  Thus,  the  above  comparison  includes  only 
design,  structural  materials,  and  labor.  It  is  still  an 
illustrative  comparison,  which  gives  an  impression  of 
potential.  Thus,  the  flexibility  which  will  be  designed  and 
built  into  the  FRBF  can  provide  for  neutron  transport, 
tritium  breeding,  and  power  density  investigations  in  a 
variety  of  blankets  at  a  relatively  low  cost.  The 
importance  of  this  can  be  seen  by  reviewing  recent 
literature  related  to  the  BCSS.  The  number  of  blanket 
design  candidates  has  been  reduced  to  only  four  due  to  the 
expenses  involved  in  testing  and  analyzing  them. 


Constructing  the  FRBF  at  Purdue  has  another  financial 
advantage;  the  laboratory  which  housed  the  FBBF  has  adequate 
space  for  installation  of  the  FRBF  and  its  electronic 
equipment.  One  of  the  shielding  walls  of  the  FBBF  could  be 
used  to  shield  the  FRBF.  and  an  opposite  (parallel)  wall  has 
only  earth  on  the  other  side.  Thus,  only  the  construction 
of  two  shielding  walls,  the  installation  of  the  neutron 
generating  and  instrumentation  system,  and  a  new  blanket  are 
required  to  have  created  a  new  fusion  reactor  blanket 
testing  facility.  This  new  facility  could  then  be  used  for 
conducting  much  needed  differential  and  integral  transport 
experiments  on  realistic  blanket  models.  Another  advantage 
of  putting  this  new  facility  in  the  same  building  as  the 
FBBF  will  be  safety. 

VII. 2.  Operational  and  Personnel  Safety  in  the  FRBF 

Four  major  safety  issues  must  be  considered  in  the  design 
of  the  FRBF  on  a  university  campus;  they  are:  1)  radiation 
leakage,  2)  tritium  leakage,  3)  experimental  hazards,  and  4) 
criticality.  Only  two  of  these  issues  are  crucial  for  the 
FRBF — radiation  and  tritium  leakage.  The  criticality  issue 
does  not  exist  here  because  no  fissile  fuel  will  be  present; 
and  the  experimental  hazards  have  been  fully  analyzed  and 
safe  techniques  developed  with  work  for  the  FBBF.  Thus, 
radiation  leakage  and  tritium  leakage  must  be  considered. 
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The  mean-f ree-path  of  14  MeV  neutrons  in  most  fusion 
reactor  materials  is  about  10  cm.19  Due  to  the  greater 
penetration  of  these  neutrons,  more  shielding  and  a  thicker 
reflector  may  be  required.  The  FBBF  was  designed  for  a 
maximum  neutron  production  rate  of  ~1011  n/s .  However,  the 
design  was  based  on  a  fission  spectrum,  not  a  fusion 
spectrum.  Theref ore , ' when  the  FRBF  final  design  is 
completed,  an  evaluation  of  the  blanket  neutronics  should 
include  analysis  of  the  leakage  flux  to  demonstrate  there  is 
no  need  for  extra  shielding  or  reflector  materials. 

However,  because  the  blanket  and  shielding  on  the  FRBF  are 
thicker  than  the  FBBF,  and  because  fusion  blankets  are 
intentionally  constructed  of  neutron  absorbing  materials 
(while  the  FBBF  contains  fissile  materials),  a  significant 
difference  in  dose  rates  should  not  exist  for  equal  source 
strengths . 

The  FBBF  has  been  thoroughly  analyzed  for  neutron  and 
gamma-ray  dose  rates  inside  the  experimental  chamber,  inside 
the  laboratory,  and  in  the  areas  above  and  around  the 
facility.42  Additionally,  dose  rate  measurements  have  been 
made  in  these  locations;  the  dose  rate  in  the  lecture  halls 
above  the  facility  is  low  enough  to  be  not  detectable.  To 
estimate  the  dose  rates  from  the  proposed  FRBF,  only  a 
comparison  of  neutron  flux  and  spectrum  at  the  outer  edge  of 
the  reflectors  is  required.  The  ratios  of  dose  rates  will 
be  approximately  equal  to  the  ratios  of  the  leakage  fluxes. 
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Special  consideration  should  be  given  to  shielding  the  gap 
provided  for  the  neutron  source  d-beam  system.  The  other 
safety  issue,  tritium  leakage,  is  more  difficult  to 
estimate . 

Tritium  will  be  contained  in  a  titanium-tritide  (TiT2) 
coating  on  a  copper  target,  in  the  center  of  a  stainless 
steel  vacuum  tube.  Tritium  atoms  migrate  rather  easily 
through  most  substances,  and  tritium  atoms  can  be  freed  from 
the  TiT  matrix.  High-energy  deuterons,  recoil  alphas  (from 
the  d-T  reaction),  high-energy  neutrons,  and  thermal  effects 
all  dislodge  the  tritium  atoms  from  the  matrix.39  This 
tritium  can  be  easily  removed  by  the  vacuum  pump  which  is 
required  for  maintaining  a  vacuum  in  the  d-beam  tube.  There 
are  many  techniques  for  pumping  and  trapping  the  hydrogen 
isotopes,  such  as  sputtei — ion  pumps,  cryogenic  pumps, 
uranium  getter  pumps,  and  catalytic  conversion  methods. 
Preventing  excess  leakage  of  tritium  should  not  be 
technically  difficult.  However,  some  tritium  will  escape 
from  the  system;  this  leakage  should  be  estimated,  at  least 
to  put  a  cap  on  the  possible  magnitude  of  concentration,  and 
to  compare  this  to  standards. 

Another  potential  hazard  will  exist  when  changing  the 
tritiated  targets.  This  will  require  opening  the  vacuum 
system;  with  the  possible  exposure  of  personnel  to  highly 
contaminated  equipment.  The  target  region  should  be 
constructed  so  that  it  can  be  removed  as  a  sealed  unit  for 


disposal,  or  for  chancing  the  target  in  a  glovebox  or  other 
safe  facility. 

Another  path  for  tritium  loss  is  from  the  breeding 
blankets.  They  will  produce  tritium  at  a  rate  comparible  to 
the  neutron  production  rate  (1011  n/sec).  However,  this 
tritium  will  most  likely  be  contained  in  sealed  systems,  and 
some  will  be  recovered  for  measurement.  When  these 
containers  are  opened,  the  possibility  exists  for  larger 
quantities  to  escape  in  a  short  period  of  time.  Safe 
experimental  and  operating  procedures  for  handling  these 
tritiated  samples  will  have  to  be  developed  and  followed. 

Since  all  other  safety  concerns  have  been  well  addressed 
during  the  design,  construction,  and  experimenting  with  the 
FBBF;  and  since  the  radiation  and  tritium  leakage  concerns 
have  been  evaluated  above;  it  appears  the  FRBF  can  be 
operated  safely  in  the  basement  of  the  Physics  Building  at 
Purdue  University.  The  next  item  which  needs  addressed  is 
the  predicted  precision  of  measurements  in  the  FRBF. 

VII. 3.  Measurement  Techniques  and  Accuracy 

Most  measurement  methods  and  equipment  will  be  the  same 
as  those  used  in  the  FBBF  (see  Chapter  III).  Much 
experience  has  been  acquired,  techniques  have  been 
thoroughly  developed  and  tested,  and  confidence  in  results 
has  been  built.  The  statistical  precision  of  the 
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measurements  depends  upon  three  parameters;  measurement 
time,  neutron  flux,  and  precision  of  the  neutron  source 
prediction  (because  results  will  be  expressed  as  an  absolute 
ratio  of  measurement/source  neutron).  The  design  of  the 
line  source  has  a  distinct  benefit  here,  it  permits  long 
measurement  times,  on  the  order  of  weeks  or  months,  without 
changes  in  source  strength,  and  with  relatively  high  neutron 
flux.  Whereas  most  neutron  generator  targets  have 
relatively  short  half-lives  due  to  high  local  d-beam  current 
and  associated  high  temperatures;  the  stretched  design  of 
the  FRBF  source  will  reduce  the  local  current  and  heating 
rate  while  allowing  for  cooling  and  recombination  of  tritium 
with  the  TiT  matrix.  The  local  heating  rate  can  be  reduced 
by  a  factor  of  100  simply  because  the  target  is  100  cm  long, 
and  normal  neutron  generators  have  beams  of  about  1  cm2. 

This  design  will  also  permit  removal  of  heat  without 
excessive  cooling  structure  or  coolant,  thus  reducing  the 
perturbation  of  the  neutron-source  spectrum.  This  should 
allow  a  target  lifetime  about  two  to  three  orders-of- 
magnitude  greater  than  normal  beam-target  systems.  By 
monitoring  the  recoil  alphas  and  automatically  adjusting  the 
d— beam  current  to  maintain  a  constant  neutron  production 
rate,  measurements  can  be  made  over  days,  weeks,  or  even 
months.  Thus,  total  fluence  can  be  increased  easily  to 
reduce  counting  errors  below  one  percent  (one  standard 
deviation).  Then,  the  measurement  results  are  mostly 
influenced  by  efficiencies  and  source-strength 


determination.  The  alpha-recoil  method  has  been  developed 
for  determining  the  strength  of  neutron  generator  sources. 

By  using  a  bank  of  these  detectors  mounted  vertically  in 
line  with  the  blanket  axis,  and  shielded  by  the  concrete 
wall  of  the  blockhouse  and  by  lead  collimators,  neutron 
source  strength  and  axial  shape  may  be  precisely  determined. 
The  overall  measurement  statistics  would  then  be  on  the 
order  of  one  to  two  percent,  which  is  much  less  than 
deviations  between  predictions  and  measurements  in  previous 
experiments . 

With  a  source  strength  of  1011  n/s ,  exposure  of  samples 
would  be  about  10*  14  MeV  n/s  (integral  flux  is  about  4 
times  this).  ~10*  n/s  in  the  inner  blanket,  and  ~104  n/s  in 
the  outer  blanket  or  reflector.  These  values  are  based  on 
the  flux  results  from  blanket  HeC.  These  relatively  high 
values  of  neutron  flux,  combined  with  a  large  volume 
available  for  experiments  (~60  cm  axially.  ~60  cm  radially, 
and  00-90°  azimutha 1 l y) ,  allows  an  increased  capability  and 
flexibility  for  measurements.  Then,  the  use  of  multiple 
measurement  techniques  and  long  measurement  times  in  a  large 
experimental  region  could  yield  the  same  results  for  the 
FRBF  that  they  have  for  the  FBBF ;  identification  of 
deficiencies  in  methods  used  to  predict  neutron  and  gamma- 
ray  transport  in  blankets,  and  identification  of  the  causes 
of  those  deficiencies  and  elimination  of  those  causes.  A 
few  other  advantages  of  constructing  the  FRBF  at  Purdue 


University  should  be  mentioned  in  this  report. 


VII. 4.  Other  Issues 

In  addition  to  the  experience  with  experimental 
techniques  previously  discussed,  there  are  other  advantages 
at  Purdue.  One  of  the  major  advantages  of  installing  the 
FRBF  at  Purdue  is  its  resources — the  people  and  the 
facilities.  Much  experience  has  been  gained  among  faculty, 
staff,  and  students  at  Purdue  while  designing,  constructing, 
and  experimenting  in  the  FBBF.  Benefits  can  be  derived  from 
this  experience,  which  includes: 

•  Theory  of  fast-neutron  reactions, 

•  Theory  of  cross  section  generation  for  transport 
predictions , 

•  Development  and  use  of  neutron  and  gamma  measurement 
techniques , 

•  Development  of  a  vectorized  integral  transport 
code , 4  4 

•  Design,  fabrication,  acquisition,  and  installation 
of  components  and  systems  in  the  FBBF,  and 

•  Analysis  and  comparison  of  fast  reactor  neutronics 
to  measurements. 

Another  benefit  is  the  physical  resources,  for  instance,  the 
Cyber  205  vector  supercomputer  is  available  for  all  faculty 
and  students,  and  for  all  research  projects.  This 
experience  and  expertise  could  contribute  significant 
improvements  to  the  technology  required  for  the  future  of 
fusion  power. 
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VIII.  CONCLUSION 


The  FRBF  is  a  conceptual  design  of  a  new  fusion  reactor 
blanket  simulation  facility  unlike  those  used  for  past 
research.  Deficiencies  in  the  designs  of  past  fusion 
reactor  blanket  neutronics  experiments  have  largely  been 
eliminated  in  the  design  of  the  FRBF.  The  neutronics 
evaluation  conducted  for  this  research  indicates  that  the 
energy-dependence  of  flux  can  be  separated  from  the  axial 
dimension  because  of  the  use  of  an  axially  distributed, 
chopped-cosine  source  with  a  cylindrical  geometry. 

Neutronics  predictions  will  produce  errors  due  to  design  of 
less  than  five  percent,  if  based  on  1-dimensional 
predictions.  Actual  analyses  will  include  2-dimensional 
treatments,  which  include  geometric  asymmetry  and  source 
anisotropy.  Accurate  treatments  of  neutronics  in  the  frbf 
will  not  require  3-dimensional  analyses,  which  are 
necessarily  done  in  course  space  and  energy.  Then,  the 
inaccuracies  of  transport  predictions  in  this  facility  (C/E 
values)  will  contribute  only  a  few  percent  to  the  deviations 
to  be  evaluated.  The  errors  in  predicting  axial  leakage 
will  contribute  the  least  to  the  overall  deviations,  due  to 
the  chopped-cosine  source  profile.  Additionally,  some  of 


these  deviations  can  be  reduced  by  reducing  the  volume  of 
void  in  the  blanket  and  by  using  the  minimum  wall  thickness 
for  vacuum  enclosures.  A  capability  to  perform  measurements 
for  long  periods  and  to  use  multiple  measurement  techniques 
improves  the  quality  of  data  which  can  be  obtained.  Thus, 
the  FRBF  can  produce  C/E  values  which  are  more  accurate  than 
past  comparisons  between  experiment  and  prediction. 

Following  is  a  summary  of  this  report. 

VIII. 1.  Summary 

The  results  of  the  few  fusion  reactor  blanket  simulation 
experiments  performed  to  date  have  shown  that  improvements 
are  needed  in  neutron  transport  codes  and  data  files,  in 
measurement  techniques,  but  most  critically  in  experiment 
design.  The  conceptual  design  of  a  new  simulation  facility, 
the  FRBF,  was  developed.  This  facility  has  an  axiailly 
distributed  neutron  source  with  a  chopped-cos ine  shape.  To 
determine  if  azimuthal  asymmetry  or  neutron  source 
anisotropy  would  severely  affect  predictions  throughout  the 
facility,  neutronics  in  two  blanket  compositions  were 
studied.  The  results  of  the  neutronics  analyses  indicated 
that  the  design  features  of  the  facility  and  its  source 
would  not  contribute  in  any  significant  way  to  inaccuracies 
in  the  prediction  of  results.  Thus,  any  deviations  between 
calculated  and  experimental  values  (C/E)  can  be  interpreted 
in  terms  of  inaccuracies  of  cross  sections  (such  as 


inelastic  scattering)  and  deficiencies  of  computational 
methods  (such  as  group  constant  generation  and  neutron 
transport  codes).  Arguments  were  then  presented  to 
demonstrate  the  advantages  of  constructing  an  FRBF-type 
facility  at  Purdue  University.  Its  cost  would  be  low 
because  it  could  be  built  adjacent  to  the  FBBF.  The  FRBF 
could  be  used  to  produce  precise  measurements  of  neutron 
fluxes  in  a  variety  of  fusion  reactor  blanket  materials  and/ 
or  designs,  at  a  relatively  low  cost.  These  measurements 
could  then  be  compared  to  design  predictions,  with 
confidence  that  the  facility  would  not  contribute 
substantially  to  deviations.  Thus,  the  FRBF  could  be  used 
to  validate  neutronics  methods  and  data  libraries,  and  could 
contribute  significantly  to  capabilities  to  predict  tritium 
breeding  rates  and  power  densities  in  the  first  wall, 
blanket,  shielding,  magnets,  and  other  components  of 
proposed  fusion  reactor  designs.  Studies  conducted  in  the 
FRBF  could  greatly  enhance  the  expanding  body  of 
methodologies  for  the  future  of  fusion  technology.  The 
following  recommendations  should  be  considered  during  the 
engineering  design  of  this  facility. 

VIII. 2.  Recommendations 

Further  studies  are  required  before  construction  of  the 
FRBF.  These  include  a  determination  of  how  to  construct  the 
blanket,  analysis  of  designs  of  targets  for  the  d-beam, 


analysis  of  designs  of  accelerator  and  beam-sweeping 
systems,  design  of  a  vacuum  system  for  the  neutron 
generator,  and  a  thorough  analysis  of  shielding 
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requirements.  Some  suggestions  are: 

1)  Economical  design  of  the  blanket  should  stress 
maximum  flexibility.  It  should  provide  for  removal  of  large 
or  small,  axial,  radial,  or  azimuthal  segments.  Excessive 
radia 1 1 y— directed  channels  should  be  avoided,  as  these  will 
tend  to  promote  neutron  streaming. 

2)  The  target  support  should  have  minimum  structure  and 
coolant.  One  method  of  construction  might  be  to  attach  the 
TiT  foil  to  a  thin-walled  copper  tube  in  the  center  of  the 
blanket.  This  could  then  be  cooled  by  flowing  air  or 
another  gas  through  it.  If  the  tube  is  aligned  such  that 
the  center  of  the  target  is  at  the  center  of  the  blanket, 
with  the  wall  of  the  tube  curving  away  from  the  beam,  the 
support  will  contribute  little  disturbance  of  the  azimuthal 
flux  profile  around  the  60°-120°  region.  Also,  safety 
should  be  a  major  consideration  in  the  design  of  the  neutron 
generating  system,  see  Chapter  VII. 2  for  a  discussion. 

3)  The  neutron  generator  and  associated  electronics 
(source-strength  detectors,  magnets,  etc.)  should  be  located 
outside  the  shielding  wall  of  the  facility.  This  will 
protect  them  from  excessive  neutron  flux  from  the  blanket, 
will  reduce  background  counts  in  the  detectors,  and  will 
reduce  asymmetries  in  the  FRBF  room. 

4)  A  zone  in  the  center  of  the  blanket  for  scattering 


w. 
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the  source  neutrons  should  be  investigated.  This  could 
promote  smoothing  of  the  source  anisotropy,  reduce  the 
energy  slightly,  to  approximate  the  plasma  source,  and 
introduce  some  angular  variation  of  the  flux  in  the  source 
region.  The  material  in  this  zone  must  have  an  elastic 
cross  section  which  is  orders-of-magni tude  greater  than  the 
inelastic  cross  section;  otherwise,  an  unrealistic  boundary 
source  will  be  imposed  on  the  first  wall,  and  the  scatterer 
will  attenuate  high-energy  neutrons  instead  of  just 
scattering  them. 

5)  The  flux  variations  found  in  the  low-energy  groups 
in  the  second  model  (Blanket  LLC ,  LiPb— cooled)  should  be 
investigated.  The  source  of  the  deviations  should  be 
determined,  and  a  method  for  ameliorating  the  effect  should 
be  found. 
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♦Read  2.749e-2  as  2.749- 10-t0  atom/barn- 


le  A. 3 


Fifty-group  Structure  (ORNL) . 


Group 

Energy 

(MeV) 

Lethargy 

Spectrum 

FRBF 

Spectrum 

1 

1 . 568e+7* 

-5 . OOe— 1 

3. 124e— 1 

1 . 069e-2 

2 

1 .492e+7 

-4. 00e-l 

2. 135e+0 

1 . 563e+0 

3 

1 . 455e+7 

-3 . 75e-l 

6 . 420e+0 

2 . 763e+0 

4 

1 . 419e+7 

-3 . 50e— 1 

8. 325e+0 

2 . 665e+0 

5 

1 . 384e+7 

-3. 25e-l 

4 . 787e-l 

2. 476e+0 

6 

1 . 350e+7 

-3. 00e-l 

1 .409e+0 

5 . 220e-l 

7 

1 . 252e+7 

-2 . 25e-l 

9. 205e-4 

0.000 

8 

1 . 221 e+7 

-2 . OOe— 1 

3 . 682e-3 

0.000 

9 

1 . 105e+7 

-1.00e-l 

3 . 682e-3 

0.000 

10 

1 . OOOe+7 

0 . 00e+0 

4. 882e— 3 

0.000 

11 

9. 048e+6 

1 . 00e-l 

1 . 234e-3 

0.000 

12 

8. 187e+6 

2 . OOe— 1 

1 . 234e-2 

0.000 

13 

7 . 408e+6 

3 . OOe— 1 

1 .801e-2 

0.000 

14 

6 . 703e+6 

4 . OOe— 1 

2 . 500e-2 

0.000 

15 

6 . 065e+6 

5. OOe-1 

3 . 315e— 2 

0.000 

16 

5 . 488e+6 

6 . OOe-1 

4. 221e— 2 

0.000 

17 

4. 966e+6 

7 . OOe— 1 

5. 129e-2 

0.000 

18 

4. 493e+6 

8. OOe-1 

6. 143e— 2 

0.000 

19 

4. 066e+6 

9. OOe-1 

7 . 069e— 2 

0.000 

20 

3 . 679e+6 

1 ,00e+0 

7 . 912e-2 

0.000 

21 

3. 329e+6 

1 . 10e+0 

8. 368e— 2 

0.000 

22 

3 . 012e+6 

1 .20e+0 

9. 220e-2 

0.000 

23 

2. 725e+6 

1 .30e+0 

9 . 643e-2 

0.000 

24 

2 . 466e+6 

1 . 40e+0 

9 . 900e— 2 

0.000 

25 

2. 231e+6 

1 . 50e+0 

9. 994e-2 

0.000 

26 

2. 019e+6 

1 . 60e+0 

1 . OOOe-1 

0.000 

27 

1 . 827e+6 

1 .70e+0 

1 . OOOe-1 

0.000 

28 

1 . 653e+6 

1 . 80e+0 

1. OOOe-1 

0.000 

29 

1 . 496e+6 

1 . 90e+0 

1. OOOe-1 

0.000 

30 

1 . 353e+6 

2. 00e+0 

2. OOOe-1 

0.000 

31 

1 . 108e+6 

2. 20e+0 

4. OOOe-1 

0.000 

32 

7 . 427e+5 

2. 60e+0 

3. OOOe-1 

0.000 

33 

5 . 502e+5 

2 . 90e+0 

3. OOOe-1 

0.000 

34 

4. 076e+5 

3. 20e+0 

3. OOOe-1 

0.000 

35 

3. 020e+5 

3. 50e+0 

2 . 500e-2 

0.000 

*Read  1 

. 368e+7  as  1.568- 10l°. 

Table  A. 3,  Continued 


Group 

Energy 

(Me/) 

Lethargy 

Spectrum 

FRBF 

Spectrum 

36 

2 . 945e+5 

3 . 53e+0 

7 . 500e-2 

0.000 

37 

2 . 732e+5 

3. 60e+0 

1.000e-l 

0.000 

38 

2. 472e+5 

3 . 70e+0 

1 . 000e-l 

0.000 

39 

2 . 237e+5 

3. 80e+0 

1.000e-l 

0.000 

40 

2. 042e+5 

3. 90e+0 

1 . 000e-l 

0.000 

41 

1 .832e+5 

4. 00e+0 

1.000e-l 

0.000 

42 

1 . 657e+5 

4. 10e+0 

1 . 000e-l 

0.000 

43 

1 . 500e+5 

4. 20e+0 

8. 000e-l 

0.000 

44 

6 . 738e+4 

5 . 00e+0 

1 . 000e-l 

0.000 

45 

2. 479e+4 

6. 00e+0 

2. 000e+0 

0.000 

46 

3. 355e+3 

8. 00e+0 

2. 000e+0 

0.000 

47 

4 . 540e+2 

1 . 00e+l 

2. 000e+0 

0.000 

48 

6. 142e+l 

1 . 20e+l 

2. 000e+0 

0.000 

49 

8. 315e+0 

1 . 40e+l 

3 . 581e+3 

0.000 

50 

1 . 000e-l 

1 .84e+l 

3. 161e+4 

0.000 

Table  A. 4.  Forty-one-group  Structure  (Purdue). 


Group 

Energy 

(MeV) 

Lethargy 

Spectrum 

FRBF 

Spectrum 

1 

1 . 568e+7* 

-5. 00e-l 

2. 447e+0 

1 . 561e+0 

2 

1 .455e+7 

-3 . 75e-l 

6 . 420e+0 

2 . 763e+0 

3 

1 . 419e+7 

-3 . 50e-l 

8 . 325e+0 

2. 665e+0 

4 

1 . 384e+7 

-3 . 25e-l 

2. 787e— 1 

2. 476e+0 

5 

1 . 350e+7 

-3. 00e-l 

1 . 409e+0 

5. 220e-l 

6 

1 . 252e+7 

-2 . 25e-l 

8. 285e— 3 

0.000 

7 

1 . 000e+7 

0. OOe+O 

6 . 822e-2 

0.000 

8 

6. 065e+6 

5. 00e-l 

2. 593e— 1 

0.000 

9 

3 . 679e+6 

1 . 00e+0 

2. 531e-l 

0.000 

10 

2 . 231 e+6 

1 . 50e+0 

5. 000e-l 

0.000 

11 

1 . 353e+6 

2. 00e+0 

5. 000e-l 

0.000 

12 

8. 208e+5 

2. 50e+0 

6. 000e-l 

0.000 

13 

4. 505e+5 

3. 10e+0 

4. 000e-l 

0.000 

14 

3. 020e+5 

3. 50e+0 

5. 000e-l 

0.000 

IS 

1 .832e+S 

4. 00e+0 

5. 000e-l 

0.000 

16 

1 . llle+5 

4 . 50e+0 

5. 000e-l 

0.000 

17 

6 . 738e+4 

5. 00e+0 

5. 000e-l 

0.000 

18 

4. 087e+4 

5 . 50e+0 

5. 000e-l 

0.000 

19 

2 . 479e+4 

6. 00e+0 

5. 000e-l 

0.000 

20 

1 .503e+4 

6 . 50e+G 

5. 000e-l 

0.000 

21 

9. 1 19e+3 

7 . 00e+0 

5 . 000e-l 

0.000 

22 

5 . 531 e+3 

7 . 50e+0 

5 . OOOe— 1 

0.000 

23 

3. 355e+3 

8. 00e+0 

5 . 000e-l 

0.000 

24 

2. 035e+3 

8. 50e+0 

5 . OOOe— 1 

0.000 

25 

1 .234e+3 

9. 00e+0 

5 . OOOe— 1 

0.000 

26 

7. 485e+2 

9 . 50e+l 

5 . OOOe— 1 

0.000 

27 

4 . 540e+2 

1 . 00e+l 

5. 000e-l 

0.000 

28 

2 . 754e+2 

1 . 05e+l 

5 . OOOe— 1 

0.000 

29 

1 . 670e+2 

1 . 10e+2 

5. 000e-l 

0.000 

30 

1 . 013e+l 

1 . 15e+l 

5 . 000e-l 

0.000 

31 

6. 144e+l 

1 . 20e+l 

5 . OOOe— 1 

0.000 

32 

3 . 727e+l 

1 . 25e+l 

5. 000e-l 

0.000 

33 

2. 260e+l 

1 . 30e+l 

5 . 000e-l 

0.000 

34 

1 . 371 e+1 

1 . 35e+l 

5. 000e-l 

0.000 

35 

8. 315e+0 

1 . 40e+l 

5. 000e-l 

0.000 

36 

5. 043e+0 

1 . 45e+l 

5. 000e-l 

0.000 

37 

3 . 059e+0 

1 .50e+l 

5 . 000e-l 

0.000 

38 

1 . 855e+0 

1 . 55e+l 

5 . 000e-l 

0 . 000 

39 

1 . 125e+0 

1 . 60e+l 

5 . OOOe- 1 

0.000 

40 

6. 826e-l 

1 . 65e+l 

5. 000e-l 

0.000 

41 

4 . 140e+0 

1 . 70e+l 

3 . 519e+4 

0.000 

*Read  1 

. 568e+7  as  1.568- 10‘°. 

V"V  .» 
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Appendix  B.  Angle-energy  Correlation  Code 

The  following  program  reads  dE/dx  data  for  Ti  and  T  and 
neutron-production  cross  sections  for  the  T(d,n)4He 
reaction.  It  then  computes  the  combined  energy  loss  rates 
for  TiT  and  does  a  numerical  integration  of  the  slowing  down 
equation.  It  uses  an  angle-energy  relation  based  upon 
reaction  kinematics  to  produce  a  table  of  energy  group  vs 
emission  angle  of  d—T— generated  neutrons.  It  also  predicts 
the  neutron  yield  per  ampere  of  d-beam  current.  This  code 
is  derived  from  ORNL/TM-9251 ;  however,  the  integration 
method  has  been  changed  to  better  approximate  the  results  in 
the  forward  and  backward  directions. 
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